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Summary. - We give in this lecture an introduction to the physics of two- 
dimensional (2d) Bose gases. We first discuss the properties of uniform, infinite 
2d Bose fluids at non-zero-temperature T. We explain why thermal fluctuations 
are strong enough to destroy the fully ordered state associated with Bose-Einstein 
condensation, but are not strong enough to suppress superfluidity in an interact- 
ing system at low T. We present the basics of the Berezinskii-Kosterlitz-Thouless 
theory, which provides the general framework for understanding 2d superfluidity. 
We then turn to experimentally relevant finite-size systems, in which the presence 
of residual "quasi-long-range" order at low temperatures leads to an interesting in- 
terplay between superfluidity and condensation. Finally we summarize the recent 
progress in theoretical understanding and experimental investigation of ultra-cold 
atomic gases confined to a quasi 2d geometry. 
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1. Introduction 

The properties of phase transitions and the types of order present in the low-tempera- 
ture states of matter are fundamentally dependent on the dimensionality of physical 
systems [1]. Generally, highly ordered states are more robust in higher dimensions, while 
thermal and quantum fluctuations, which favor disordered states, play a more important 
role in lower dimensions. 

The case of a two-dimensional (2d) Bose fluid is particularly fascinating because of 
its "marginal" behavior. In an infinite uniform 2d fluid thermal fluctuations at any non- 
zero temperature are strong enough to destroy the fully ordered state associated with 
Bose-Einstein condensation, but are not strong enough to suppress superfluidity in an 
interacting system at low, but non-zero temperatures. Further, the presence of residual 
"quasi-long-range" order at low temperatures leads to an interesting interplay between 
superfluidity and condensation in all experimentally relevant finite-size systems. 

This behavior is characteristic of a wide range of physical systems which share some 
generic properties such as dimensionality, form of interactions, and Hamiltonian sym- 
metries. These include liquid helium films [2], spin-polarized hydrogen [3], Coulomb 
gases [4], ultra-cold atomic gases [5, 6], exciton [7, 8] and polariton [9, 10] systems. More- 
over, the Berezinskii-Kosterlitz-Thouless theory which provides the general framework 
for understanding 2d superfluidity is also applicable to a range of physical phenomena 
in discrete systems, such as ordering of spins on a 2d lattice [11] and melting of 2d 
crystals [12, 13]. 

In this paper we give an introduction to the physics of 2d Bose fluids from an atomic 
physics perspective. Our goal is to summarize the recent progress in theoretical under- 
standing and experimental investigation of ultra-cold atomic gases confined to 2d geom- 
etry, and we also hope to provide a useful introduction to these systems for researchers 
working on related topics in other fields of physics. 

IT. Absence of true long-range order in 2d. - The most familiar phase transitions 
in three spatial dimensions (3d), such as freezing of water, ferromagnetic ordering in 
spin systems, or Bose-Einstein condensation (BEC), are all associated with emergence 
of true long-range order (LRO) below some non-zero critical temperature. Such order is 
embedded in a spatially uniform order parameter, e.g. magnetization in a ferromagnet 
or the macroscopic wavefunction ip describing a BEC. Further, in all the above cases 
emergence of true LRO corresponds to spontaneous breaking of some continuous sym- 
metry of the Hamiltonian. In case of crystallization (freezing) , translational symmetry is 
spontaneously broken. For Heisenberg spins on a (fixed) lattice, the Hamiltonian has a 
continuous spin rotational symmetry which is spontaneously broken in the ferromagnetic 
state. In case of BEC, the phase of tp is arbitrarily spontaneously chosen at the transi- 
tion. As we will discuss in more detail later, under certain conditions this makes a Bose 
gas formally equivalent to a system of two-component spins on a lattice, the so-called 
XY model. 

Already in 1934 Peierls pointed out that the possibility for a physical system to 
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exhibit true LRO can crucially depend on its dimensionality [14, 15]. Peierls considered 
a 2d crystal at finite temperature T, and studied the effects of thermal vibrations of 
atoms around their equilibrium positions in the lattice (i.e. phonons). He found that 
the uncertainty in the relative position of two atoms diverges with the distance between 
their equilibrium positions: 



for r = |r| ^> a, where u (r) is the atom displacement from its equilibrium position r, a 
is the lattice spacing, and (...) denotes a thermal average. In contrast, the corresponding 
result in 3d is finite, and small compared to a if T is below the melting temperature. 
The result in eq. (1) is in direct contradiction with the starting hypothesis of long-range 
crystalline order, since it implies that, based on the positions of atoms in one part of the 
system, we cannot predict with any certainty the positions of atoms at large distances. 

The Peierls result (1) is a simple example of the absence of spontaneous symmetry 
breaking at non-zero T in 2d. The absence of LRO in low-dimensional systems was later 
more generally and formally studied by Bogoliubov [16], Hohenberg [17], and Mermin 
and Wagner [18]. The general statement is that LRO is impossible in the thermodynamic 
limit at any non-zero T in all Id and 2d systems with short-ranged interactions and a 
continuous Hamiltonian symmetry. This is now most commonly known as the Mermin- 
Wagner theorem. In all such systems Hamiltonian symmetry is always restored by low- 
energy long-wavelength thermal fluctuations, the so-called Goldstone modes. In the case 
of an interacting Bose gas Goldstone modes are phonons, while in the case of the XY 
model on a lattice they are spin-waves. As a direct consequence of the functional form of 
the density of states in low dimensionality, such modes always have a diverging infrared 
contribution and destroy LRO. It is however important to stress that the absence of true 
LRO does not preclude the possibility of any phase transitions in 2d systems, just the 
symmetry breaking ones. As we will see later, a phase transition associated with the 
apparition of a topological order does take place in a uniform 2d gas. 

12. Outline of the paper. - In Section 2 we will explicitly see that Bose-Einstein 
condensation is impossible in both the ideal and a repulsively interacting infinite uniform 
2d Bose gas. The absence of true LRO in these systems is seen in the fact that the first 
order correlation function: 



where \fr(r) is the annihilation operator for a particle at position r, always tends to zero 
for r — > oo. Note that these two statements are equivalent under the Penrose-Onsager 
definition of the condensate density [19]: 



(1) 




(2) 




(3) 



n = lim gi (r) 

r— >oo 
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However, the weak logarithmic divergence in eq. (1) suggests that the destruction 
of LRO is only marginal in 2d. The consequence of this weak divergence is that an 
interacting Bose gas at low T exhibits "quasi- long-range order", corresponding to g\{r) 
which decays only algebraically with distance. 

This low T state is also superfluid, and the phase transition between the superfluid 
and the high T normal state is described by the Berezinskii-Kosterlitz-Thouless (BKT) 
theory [20, 21], which we discuss in Section 3. Further, the slow decay of gi(r) at low T 
has important consequences for condensation and symmetry breaking in the experimen- 
tally relevant finite-size systems. We address this issue in Sections 4 and 5, first for a 
finite box potential, and then for the experimentally most pertinent case of a harmoni- 
cally trapped gas. 

In Sections 6 and 7 we introduce the experimental methods of atomic physics used in 
the current studies of 2d Bose gases. In § 6 we give an overview of experimental systems 
in which (quasi-)2d Bose gases have been realized, and in § 7 we discuss the experimental 
probes of coherence in these systems. We conclude by outlining some research directions 
which are likely to be of interest in the near future in Section 8. 

2. The infinite uniform 2d Bose gas at low temperature 

In this section we discuss thermal fluctuations and the absence of Bose-Einstein con- 
densation in an infinite uniform 2d gas. In § 2T we consider the ideal gas, in which no 
phase transition occurs. In this case the first order correlation function g\(r) gradually 
changes from a gaussian function in the high-temperature, non-degenerate regime, to an 
exponentially decaying function in the degenerate regime. 

In the repulsively interacting 2d Bose gas the BKT phase transition to a superfluid 
state occurs at a non-zero critical temperature, but the conclusion that BEC transition 
does not occur remains true. We first introduce the description of interactions in an 
atomic 2d gas in § 2'2. In § 2'3 we qualitatively discuss why interactions lead to a strong 
suppression of density fluctuations in a degenerate gas, so that the low-energy long- 
wavelength excitations (phonons) in this system are almost purely phase fluctuations. 
The Bogoliubov analysis presented in § 2'4 provides a more quantitative justification for 
this conclusion and also indicates why we expect an interacting 2d gas to be superfluid 
at very low T. As we show in § 2'5, the long-wavelength phonons still lead to a vanishing 
gi(r) at r — > oo. Therefore, in accordance with the Mermin- Wagner theorem, these 
"soft" Goldstone modes still destroy the long-range order and restore the Hamiltonian 
symmetry. However, the decay of g\ at large distances is only algebraic with r at very low 
T. This low temperature, superfluid state is said to exhibit "quasi-long-range order". The 
BKT phase transition from the superfluid state with algebraic correlations to the normal 
state with exponentially decaying correlations is discussed in the following Section 3. 

2T. The ideal 2d Bose gas. - The absence of Bose-Einstein condensation in the ideal 
2d Bose gas can straightforwardly be seen by following the standard Einstein's argument 
which associates condensation with the saturation of the excited single-particle states at 
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some non-zero temperature. 

In 2d, the density of states for spinless bosons mL 2 /(2irh 2 ) is constant, where m is 
the particle mass and L — > oo is the linear size of the system. Assuming no condensation, 
the total number of particles is then: 



(4) N= mL2 ^ " 



2ttH 2 ,L eP( s -ri - 1 



where j3 = 1/(/cb? 1 ), and fi < is the chemical potential. Equivalently, the phase-space 
density D is given by: 

f°° rlr 

(5) D^n\ 2 = —— = - H l-Z) , 



where n = N/L 2 is the 2d number density, A = hj \f^wmk^T is the thermal wavelength, 
and Z = e 13 ^ is the fugacity. 

In 3d, the signature of BEC is that the analogous relationship between the phase 
space density and fugacity has no solutions for Z when the phase space density is larger 
than the critical value n 3 A 3 ~ 2.612, where n 3 is the 3d number density. Below the 
condensation temperature, chemical potential is fixed at ji = and the phase space 
density of particles in the excited states is saturated at w 2.612. However, we see that 
in 2d a valid solution: 



(6) = 1 - e 



always exists. In other words, for any non-infinite phase space density there exists a 
negative value of fi which allows normalization of the thermal distribution of particles 
in the excited states to the total number of particles in the system N. This shows that 
BEC docs not occur in the ideal infinite uniform 2d Bose gas. 

We next look at the first order correlation function gi(r), which we can write as the 
Fourier transform of the momentum space distribution function n& : 

If 00 1 h 2 k 2 

(7) Qi(r) = -. — ttt / ni. e tk r d 2 k . with nu = —r, r , et = . 

y> yn ' (2tt) 2 Jo k ' fe e /3(e fc -f*) - 1 ' 2m 

In absence of condensation gi(r) always vanishes at r — > oo. However, it still shows 
qualitatively different behavior at high and low temperature: 

• In a non-degenerate gas, eq. (6) gives Z w n\ 2 <C 1. In this regime \/j,\ ^> k^T and 
all momentum states are weakly occupied: 



(8) 



nk^Ze-^x «nA 2 e- fc2A2 / 4T «1 (Vfe). 
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In this case gi(r) is gaussian, showing only short-range correlations which decay 
on the length scale \/^/tt: 



(9) ffi(r)*ne' 



-7rr 2 /A 2 



In a degenerate gas with nX 2 > 1, from eq. (6) we get Z ps 1 and « e~ IlA <C 1, 
so that D = nX 2 ss ln(/c B r/|^|). In this regime the occupation of high energy states, 
with (3ek ~> 1, is still small and given by the Boltzmann factor: 

(10) n k M e"** = e -* 2 * 2 / ** < x j for fc 2 > 4 ^ /A 2 _ 

However, the low energy states with f3ek <C 1 are strongly occupied: 

where A: c = ^/2m|/x|//i. In this case gi(r) is bimodal. At short distances, up to 
r ~ A, correlations are still gaussian as in eq. (9). However, the Lorentzian form 
in cq. (11) corresponds to approximately ( 1 ) exponential decay of g\{r) at larger 
distances, r » A: 

(12) gi (r) « , with £ = fc" 1 w A e nx2 / 2 /V^. 

We can also estimate the partial phase space densities corresponding to the gaussian 
and the Lorentzian parts of the momentum distribution, and see that most particles 
accumulate in the Lorentzian part corresponding to low momentum states: 

\2 roo 

(13) L> G ~ / n k d 2 k « 1/e « £> 

(2tt) 2 y^iv^ 

A 2 fV^/y* 



(14) ^^/ o n fc « D 

We thus see that even though there is no phase transition in this system, the first 
order correlation function gradually changes from a gaussian with short range correlations 
in the non-degenerate regime, to an exponential in the degenerate regime. Further, 
for nX 2 > 1, the correlation length £ oc e nX I 2 grows exponentially. Therefore, while 
<?i(r) formally vanishes at r — » oo at any non-zero T, the length scale over which it 
decays is exponentially large in the deeply degenerate ideal Bose gas. This has important 



(- 1 ) More precisely the Fourier transform of the 2d Lorentzian distribution l/(k 2 + k 2 ) is defined 
for r ^ and is proportional to the Bessel function of imaginary argument Ko(k c r), whose 
asymptotic behaviour is e~ k " r /\/r. 
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consequences for the finite-sized experimental systems, since for any arbitrarily large fixed 
system size L there exists a low enough non-zero temperature T for which I > L, and 
correlations span the entire system. This issue will be addressed in Sections 4 and 5. 

2'2. Interactions in a 2d Bose gas at low T. - The interaction between two atoms 
at positions r t and r j in a 3d Bose gas at low temperature is well characterized by the 
contact potential ( 2 ): 

(15) V(r i -r j ) = 4 ^a s S^ d \r i -r j ) , 

where a s is the 3d s-wave scattering length, and mechanical stability requires repulsive 
interactions o s > 0. In 2d the two-body scattering problem is in general more complicated 
and the scattering amplitude is energy-dependent [25]. We will address this issue in 
more detail in Section 6'2. However, in all experimentally relevant situations so far, the 
analysis of interactions is simplified by the fact that while the gas is kinematically 2d, 
the interactions can be described by 3d scattering. The condition for this simplification 
is that the thickness of the sample £o is much larger than the 3d scattering length a s . In 
all current atomic experiments the ratio £o/a s is larger than 30. In this case we can to 
a very good approximation write the interaction energy as: 

(16) E int = 9 -Jn 2 {r)d 2 r, 

where g is the energy-independent interaction strength and n(r) is the local density. Note 
that here and in the following we treat the density n(r) - and later the phase 0(r) - as a 
classical function. This will notably simplify the mathematical aspects of our approach, 
while capturing all important physical consequences related for example to quasi-long 
range order and to the normal-superfluid transition. 

In 3d, the interaction strength </ 3d ) = (4ttH 2 /m)a s explicitly depends on the scatter- 
ing length. However, we see on dimensional grounds that in 2d we can write: 



(17) 

so that g is a dimensionless coupling constant. The 2d healing length, which gives the 
characteristic length scale corresponding to the interaction energy, is then given by: 

(18) £ = h/yjmgn = l/y/jjn . 



( 2 ) Strictly speaking, the delta function in 3d must be regularized by using for example the 
pseudo-potential [22]. The extension of the notion of zero-range potential to the 2d case is 
discussed in [23] (see also [24] for a discussion in terms of many-body T matrix) . 
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We can anticipate on dimensional grounds that g <~ o s /£q, corresponding to n = n^lo. 
Specifically, if a gas is harmonically confined to 2d, say in the x — y plane, we will see in 
Section 6'2 that the expression for g is: 

(19) ~9 = V8n°^, 

where in this case £o = a z is the oscillator length along the kinematically frozen direction 
z. 

We can qualitatively define the strongly interacting limit by the value of g for which 
the interaction energy of N particles, E lnt , reaches the kinetic energy E K of N non- 
interacting particles equally distributed over the lowest N single-particle states. In this 
limit, we expect the many-body ground state to be strongly correlated. We can estimate 
the interaction energy with a mean-field approximation by setting n(r) — n in eq. (16), 
which corresponds to neglecting density fluctuations: 

(20) E int = ^~gNn . 

2m 

Using the 2d density of states, mL 2 / (2ttH 2 ), the energy of the JV-th excited single-particle 
state is En = (2irh 2 /m)n, and so: 

(21) E K = InE n = —Nn . 

2 m 

The strongly interacting limit then corresponds to: 

(22) E int =E K ^g = 2Tr. 

For comparison, the value of g in the current experiments on atomic 2d Bose gases varies 
between <~ 10 _1 [26, 27] and ~ 1CP 2 [28], while in the more strongly interacting 4 He 
films [2] it is estimated to be of order 1 [29] . 

The fact that result (22) is independent of the density of the gas n is a natural 
consequence of the fact that g is dimcnsionless. This is in contrast to the 3d case, where 
the relative importance of interactions is characterized by the dimensionless parameter 
n 3 a 3 s . 

2'3. Suppression of density fluctuations and the low-energy Hamiltonian. - At strictly 
T = 0, a weakly interacting 2d Bose gas is condensed and described by a constant macro- 
scopic wavefunction (a uniform order parameter) ip = ^/ne %e , where n and 9 are classical 
fields. At any finite T, both the amplitude and the phase of ip show thermal fluctuations. 
However, repulsive interactions will always lead to a reduction of density fluctuations in a 
low temperature gas. The interaction energy is (g/2) J n 2 (r) d 2 r = (g/2)L 2 (n 2 (r)), and 
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so keeping the average density n = (n(r)) fixed, we see that minimizing the interaction 
energy is equivalent to minimizing the density fluctuations: 

(23) (An) 2 = (n 2 (r)) - n 2 = (g 2 (0) -l)n 2 , 

where g 2 (r) = (n(r)n(O)} /n 2 is the normalized second order (density-density) correlation 
function ( 3 ). In the ideal Bose gas 32(0) = 2, while if the density fluctuations are 
completely suppressed 32(0) = 1. 

We can estimate the minimum cost of density fluctuations by equating it with the 
increase in interaction energy from adding a single particle to the system: 

(24) — - = gn = —gn , 

oN m 

where we have set 32(0) = 1, which is consistent with the density fluctuations being 
significantly suppressed. Comparing this with the thermal energy k^T, we get: 

(25) f= = f D . 
y ' k B T 2tt 

This strongly suggests that at sufficiently low temperature, given by D 3> 27r/g, any 
density fluctuations must be very strongly suppressed. Numerical calculations [29] show 
that they can be significantly suppressed already for D ^> 1, with the exact extent of the 
suppression depending on the strength of interactions g. 

In the limit of strong suppression of density fluctuations, the interaction energy be- 
comes just an additive constant ((l/2)gn 2 L 2 ) in the Hamiltonian, and the kinetic energy 
arises only from the variations of 9, the phase of tp. The system is then often described 
by an effective low-energy Hamiltonian: 

(26) H e = ^-n s J (W) 2 d 2 r, 

where one heuristically replaces the total density n with the (uniform) superfluid density 
n s < n. This is physically motivated because one expects only the superfluid component 
to exhibit phase stiffness and to flow under an imposed variation of 8, with local velocity 
of the superfluid given by v s — (fi/m)V0. Also note that at T — superfluid density is 
equal to the total density, and at very low T they are similar. In essence, renormalizing n 
to the lower n s is an effective way of absorbing all the short distance physics, including any 
residual density fluctuations, and He provides a good description of long-range physics, 
at distances r » £, A. 



( 3 ) Note that following the conventions in the literature on different topics we normalize 32 so 
that it is dimensionless, while gi has units of density. 
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The effective low-energy Hamiltonian Hg is the continuous version of the Hamiltonian 
of the XY model of spins on a lattice. It can be used to derive the correct long-range 
algebraic decay of gi(r) in the low T superfluid state (see § 2'5). However, it is important 
to stress some caveats: 

(1) Even though at low temperature n s s=y n, Hg fundamentally cannot be the correct 
microscopic Hamiltonian. We can see this on very general grounds since the proper 
Hamiltonian is by definition a temperature-independent entity, while Hg depends on the 
temperature through n s . More accurately Hg represents the increase of the free energy 
of the gas if one imposes the superfluid current (H/m) V#, for instance by setting the gas 
in slow rotation (see the appendix in [6] for more details). 

(2) If and only if the density fluctuations are completely absent can the 2d Bose 
gas be formally mapped onto the XY model. This condition is essentially fulfilled for 
D ^> 2-K/g, but it is not satisfied at the BKT critical point. As discussed in Section 3, 
BKT transition occurs at a critical phase-space density D c = ln(380/.g), which for the 
experimentally relevant values of g corresponds to D c <~ 6 — 10. At that point density 
fluctuations arc significantly suppressed, but cannot be completely neglected. This is 
one of the key reasons that makes the proper microscopic theory of the BKT transition 
in the Bose gas difficult. 

In summary, Hg can be used as an effective quantum Hamiltonian for the derivation 
of some key features of long-range physics at low temperature, where n s <~ n. However 
Hg alone is not sufficient for the proper derivation of the BKT transition. 

2'4. Bogoliubov analysis. - By performing Bogoliubov analysis near T = we can see 
more explicitly why the density fluctuations arc suppressed at sufficiently low T, but the 
phase fluctuations are not, and also why it is natural to expect the low temperature state 
to be superfluid. 

The basic idea is that we start with the assumption that the T = state of a weakly 
interacting gas is described by the uniform order parameter ip = ^/ne l6 , find the excita- 
tion spectrum, and then consider the effects of the thermal occupation of the excitation 
modes. This approach may not seem justified in 2d, because in the end we find that 
thermal fluctuations destroy the order parameter, and thus invalidate our starting as- 
sumption. However we can qualitatively argue that it still works well as long as we have 
a well defined local order parameter, which is true if the order parameter is destroyed 
only at very large distances by the long-wavelength phase fluctuations. The applicability 
of the Bogoliubov approach to 2d quasi-condensates was formally justified in [30, 31]. 

Assuming contact interactions, the classical field Hamiltonian is given by 



The dynamics of ip(r,t) are governed by the Gross-Pitaevskii equation [32, 33]: 



(27) 




(28) 




Two-dimensional Bose fluids: An atomic physics perspective 



11 



which we can derive from eq. (27) by treating ip and ip* as canonical variables. 

We now define the local phase 9{r,t) from ip(r,t) = \tp(r,t)\ e ie ( r '*). Assuming that 
the density fluctuations are small, we write |^(t, t) | 2 = n(l + 2r](r,t)), with r\ <C 1 and 
J rjcPr — 0. Within this approximation we get, up to the additive constant gn 2 L 2 /2, 



(20) H = ^-n I (W(r)) 2 d 2 r + 



2m 



n(V?7(r)) 2 + 2gn 2 {rj(r))' 1 



d 2 r , 



and a set of coupled linear equations for the time evolution of 9(r, t) and r](r, t): 
(30) 



(31) 



89 

at 

drj 



2m 



v 2 e. 



It is now convenient to Fourier expand the phase and the density fields: 
(32) 9(r, t) = ]T c k (t) e ik r , r,(r, t) = ^ d k (t) e ik r , 



where k = 2ir(j x , j y )/ L, with j x ,j y integers, is a discrete variable since we consider for 
the moment a sample of finite size L 2 . We will let L — >• oo at the end of the calculation. 
The functions and r] are real, which implies c* k = c_fc and d* k = c?-fe- In addition the 
conservation of particle number J rjd 2 r = leads to d = 0. This yields the Hamiltonian 



(33) 



H = nL 2 J2 



2 1.2 



K 2 k 

2m 



\CkV + 



21,2 



K 2 k 



2m 



+ 2gn) \d k \ 



and the coupled equations of motion (for fe ^ 0): 
(34) 
(35) 



fhk 2 

Ck = -\2^ 



dk 



2m 



Ck 



For k = the equation of motion deduced from eq. (30), c = —gn/h, simply gives the 
time evolution of the global phase of the gas. 

In this formulation Ck and dk are conjugate dimensionless quadratic degrees of freedom 
corresponding to the phase and the density fluctuations, respectively. At each k we have 
a harmonic oscillator-like Hamiltonian, and from the equations of motion we can read 
off the eigenfrequencies 



(36) 



hk 2 


fhk 2 


2gn\ 


/ 2m 


,2m " 


h x) 
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which is the well known Bogoliubov result. At low k, the eigenmodes are phonons 
with uj k = ck, where c = yjgn/m. At high k we have free particle modes with uj k = 
Hk 2 1 (2m) + gn/h. The crossover between the two regimes is at k <~ l/£ = \/gn. 

From this analysis, and specifically the results in cqs. (33) and (36), we can draw 
several conclusions: 

(1) From the dispersion relation w^. we see that the excitation modes in this system 
have a non-zero minimal speed c. From the Landau criterion, we thus expect it to be 
superfluid with a critical velocity v c — c. This argument is identical to the 3d case, and 
relies just on the existence of a (reasonably) well defined local order parameter, rather 
than true LRO. We note however that the Landau criterion is a necessary, but not a 
sufficient condition for the existence of a superfluid state, because it does not address 
the question of metastability of the superfluid flow (sec e.g. [34]). 

(2) Both in the classical and in the quantum regime the harmonic oscillator at thermal 
equilibrium obeys the virial theorem, (mw 2 i 2 ) = (p 2 /m), where x and p are position and 
momentum, respectively. Reading-off the coefficients in front of |cfc| 2 and \dk\ 2 in eq. (33), 
we see that in our case this corresponds to: 



(\d k \ 2 ) = H 2 k 2 /2m 
(|c fc | 2 ) h 2 k 2 /2m + 2gn ' 



We thus explicitly see that long-wavelength phonons involve only phase fluctuations, 
since (\dk\ 2 ) <C (|cfc| 2 ) for k — > 0. On the other hand, the high-fc free particles involve 
both phase and density fluctuations in equal parts, since (|dfc| 2 ) = (|cfc| 2 ) for k — > oo. 

(3) We also explicitly see that density fluctuations are not "soft" Goldstone modes, 
because their energy cost does not vanish in the k — > limit. The effect of interactions is 
to suppress (compared to the ideal gas) the density fluctuations at length scales > £. The 
fluctuations on short length scales (< £, corresponding to k > l/£) are not suppressed 
by interactions, but those in any case do not have a diverging contribution. 

(4) On the other hand, the energy cost of phase fluctuations (phonons) vanishes for 
k — > 0. This conclusion would be the same in 3d, but the crucial difference is that in 
2d the density of states at low k leads to a diverging effect of these fluctuations and the 
destruction of true LRO, as shown in the next subsection § 2'5. 

(5) We can use this Bogoliubov analysis to provide an estimate of the relative den- 
sity fluctuations An 2 /n 2 = 4(r? 2 ). In principle the thermal equilibrium distribution of 
the Bogoliubov modes is given by the Bose-Einstein distribution. For simplicity we ap- 
proximate it in the following way; we suppose that the modes with frequency lower 
than ksT/h have an average energy k-QT/2 (classical equipartition theorem) and that 
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the higher energy modes have a negligible population. Therefore we take ( 4 ) 

(38) nL> ( + 2 g n) (\d k \ 2 ) = ( k f' 2 > * < ^ ' 

y J \2m y J Xl 17 \0, if Hu) k >k B T. 

In addition we assume that k B T is notably higher than the interaction energy gn, which 
is the case in most experiments realized with cold atoms so far. The cutoff huk = k B T 
then lies in the free-particle part of the Bogoliubov spectrum, at the wavevector kx ~ 
\fmk~ B T jh. We can now estimate the relative density fluctuations 

An 2 .^„, |2 \ L 2 f 4 k B T/2 



and we find 

(40) ~ ^ 

For realistic values of the ratio k B Tjgn (i.e. not exponentially large), we recover here 
the previously announced result that density fluctuations are notably suppressed when 
n\ 2 > 1. 

Based on these arguments, with some more quantitative justification we arrive at the 
same conclusions as before. First, for effectively describing the low temperature state of 
the gas, the most important part of the Hamiltonian in eq. (33) is the term corresponding 
to the phase fluctuations. Second, if we keep only that term, we also have to keep in mind 
that we are neglecting short-distance physics and effectively introducing a momentum 
cutoff at fc max = l/£. As before we can heuristically incorporate short-distance physics 
by replacing n — > n s to obtain: 

(41) H*n s L 2 ]T ^\^\ 2 = ^j(ye{r)Y<P ri 

k 

which coincides with the Hamiltonian Hg given in eq. (26). 

2'5. Algebraic decay of correlations. - To derive the low-T behavior of the one-body 
correlation function gi(r) = (tp* (r)ip(O)) at large distances, r 3> £, A, we start with the 
Hamiltonian Hg, and the wave function with no density fluctuations ip(r) = ^/nle %e ^ r \ 
Note that this normalization of ip leads in r — to the incorrect value g\ (0) = n s , whereas 
it should be ffi(0) = n > n s . However, as discussed earlier, density fluctuations at short 



( 4 ) We use here that c fe and d h are complex amplitudes, and that the modes k and — k are 
correlated, so that oLfc = d* k . 
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distances lead to a more complicated decay of gi , and at large r replacing n — > n s is the 
appropriate normalization. 

The long distance (r » £) behavior of g\ (r) essentially depends on the population of 
phonon modes with wave vector k <C which coincides with the momentum cutoff 
introduced in Hg. The occupation of eigcnmodes is in general given by the standard Bose 
result (exp((3fkdk) — Assuming as in eq. (40) that k^T > gn s , the phonon modes are 
in the regime hujk <C k-gT and the occupation number simplifies into k-QT/(fiWk) which 
leads to (classical equipartition theorem) 

, . 7-9^ ^ /I i2\ ^B-^ 1 

(42) phonon modes: n s L {\c k \ ) = . 

2m 2 

Introducing the real and imaginary part of the Fourier coefficients c& = c' k + ic' k , we have: 

(43) 04I 2 } = <KI 2 ) = • 

We recall that c' k and c' k are independently fluctuating variables ((c' k c' k ) = 0) and that 
the mode k and — k are correlated because 9 is real: c' fe = c'_ k and = — c"_ k . 
We want to calculate: 

(44) 9l (r) = (</>* (r)^(0)) = n s (e«'M-'«>») , 
where: 

(45) 9{r) - 9(0) = ^ c' k (cos(fc ■ r) - 1) - sin(fe • r) . 



For each independent Gaussian variable u, (e itl ) = e 2 (" ). Using eq. (43), and trans- 
forming the discrete sum over A; into L 2 /(Air 2 ) J d 2 k we obtain: 

1 y 1 - cos(fc • r) 2 



(46) ft(r)=na exp^-^ ? y p rffcj . 

The integral in the exponent has significant contributions only from modes > 1/r so 
that 1 — cos(fc • r) <~ 1. Since we restrict our analysis to r ^> A, this is not inconsistent 
with the classical field approximation which requires k < 1/A. The upper limit of the 
integral is set by the short-distance cutoff fc max = We thus expect the integral to 
be ~ ln(r/£). More formally we can note that V 2 J* (1 - cos(fc • r))k~ 2 d 2 k = (2Tr) 2 S(r), 
from which we infer: 



(47) 



f 1 - cos(fc • r) , 2 , /r\ 
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This leads to: 



( 



) 



l/(n 3 A 2 ) 



(48) 



5i (r) = n s 



r 



Note that depending on the relative size of A and £ we can also set the upper limit of the 
integral to 1/A, but this difference in the short-distance cutoff does not affect the main 
conclusion about the power law decay of correlations at large distances. 

To summarize, we have shown that in an interacting 2d Bose gas at low T, the 
first order correlation function gi(r) decays algebraically with r at large distances. The 
conclusion that g\ (r) vanishes for r — > oo is consistent with the Mermin- Wagner theorem, 
i.e. the absence of BEC and true LRO at any non-zero T. However the decay of gi(r) is 
very slow and the system exhibits a "quasi-long-range order" . Further, as discussed in 
the following Section 3, the exponent l/(n s A 2 ) is never larger than 1/4 in the superfluid 
state, making the decay of g\{r) extremely slow. The superfluid state with suppressed 
density fluctuations can be viewed as a superfluid "quasi-condensate" , i.e. a condensate 
with a fluctuating phase [35, 36, 37]. 

3. The Berezinskii Kosterlitz Thouless (BKT) transition in a 2d Bose gas 

Our analysis so far does not explain how the phase transition from the low temperature 
superfluid state to the high temperature normal state takes place. This transition is 
unusual because it does not involve any spontaneous symmetry breaking in the superfluid 
state, and in the usual classification of classical phase transition is termed "infinite order" , 
suggesting that most thermodynamic quantities (except for example superfluid density) 
vary smoothly at the transition. There is no true LRO on either side of the transition, 
but the functional form of the decay of gi(r) changes from algebraic in the superfluid 
state (corresponding to quasi-LRO) to exponential in the normal state (corresponding 
to no LRO). 

The microscopic theory of the 2d superfluid transition was developed by Berezin- 
skii [20] and Kosterlitz and Thouless [21] (see [4] for a more recent review). The transition 
takes place in the degenerate regime, where the density fluctuations in an interacting gas 
are significantly suppressed. We therefore expect that the transition can still be at least 
qualitatively explained by considering only phase fluctuations. However, a sudden tran- 
sition with a well defined critical point cannot be explained by considering only phonons, 
since we have seen in eq. (48) that, while they destroy true LRO at any non-zero T, their 
effect grows smoothly with temperature. 

3T. The role of vortices and topological order. - The key conceptual ingredient of the 
BKT theory is that in addition to phonons described by the Hamiltonian (26), another 
natural source of phase fluctuations are vortices, points at which the superfluid density 
vanishes, and around which the phase 9 varies by a multiple of 2tt. For our purposes we 
can consider only "singly-charged" vortices with phase winding ±27r, which are energet- 
ically stable. Around an isolated single vortex, centered at the origin, the velocity field 
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Fig. 1. - The BKT mechanism at the origin of the superfluid transition. Below the transition 
temperature (left figure), vortices exist only in the form of bound pairs formed with vortices of 
opposite circulation. When approaching the transition point the density of pairs grows and the 
average size of a pair diverges. Just above the transition point (right figure), a plasma of free 
vortices is formed and the superfluid density vanishes. 



hS78/m varies as h/(mr), corresponding to angular momentum h per particle. The two 
signs of the vortex charge correspond to the two senses of rotation around the vortex. 
The size of the vortex core (hole in the superfluid density) is set by the healing length 
£, so their presence is not inconsistent with the picture that the density fluctuations are 
suppressed at length scales r 3> £. In fact, it makes sense to speak of well defined indi- 
vidual vortices only if away from the vortex cores the density fluctuations are suppressed 
on the length scale £; otherwise we simply have a fully fluctuating thermal gas. 

As we will illustrate below, once one considers vortices as another source of phase 
fluctuations, one can explain the microscopic mechanism behind the superfluid-to-normal 
phase transition. Below a well defined critical temperature Tbktj vortices can exist only 
in the form of bound ("dipole") pairs of vortices with opposite circulations ±27r. Since 
they do not have any net charge, such pairs do not create any net circulation along closed 
contours larger than the pair size, which for a tightly bound pair is also of order £ ( 5 ). 
Such pairs therefore have only a short-range effect on the phase 9 and the velocity field, 
and do not have a large effect on the behavior of g±(r) at large distances. They fall under 
"short-range physics" of the system and together with the residual density fluctuations 
they lead to some renormalization of n s , but do not qualitatively alter the phenomenol- 
ogy of the long-range physics discussed in Section 2. On the other hand, above Tbkt 
unbinding of vortex pairs and proliferation of free vortices becomes energetically favor- 
able. Free vortices then form a disordered gas of phase defects and completely "scramble" 
the phase 9 (see figure 1). This destroys the quasi-LRO and suppresses superfluidity. At 
even higher temperature where density fluctuations are strong, the notion of individual 
vortices becomes physically irrelevant. 

In hindsight, we can associate superfluidity with presence of a "topological order". 
Long- wavelength phase fluctuations (phonons) destroy true LRO, but do not alter the 



( ) For a dipole field [V0| ~ 1/r 2 , so the circulation § V6* • dr vanishes for large contours. 
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topology of the system. In other words, phonons lead to smooth local variations of the 
field V which can be eliminated (or "ironed out") by continuous deformations. The 
same argument holds for bound vortex pairs which can be annihilated. Therefore, the 
superfluid quasi-condensate with no free vortices is topologically identical to the BEC 
with true LRO. On the other hand, an isolated free vortex cannot be unwound and 
eliminated from the system by continuous deformations of if); it affects the phase 9 non- 
locally, at arbitrary large distances. The annihilation argument also does not work for a 
plasma of free vortices, because if we consider a closed contour of arbitrary large size it 
will in general not contain equal number of vortices with opposite charges. Therefore at 
any length scale the free vortex plasma is topologically different from an ordered BEC. 
Although we could not have necessarily anticipated this, we can deduce that topological 
order, rather than true LRO, is a sufficient condition for superfluidity in an interacting 
2d Bose gas. 

3'2. A simple physical picture. - The full thermodynamic description of the 2d gas, 
including the role of vortices, is a difficuly task. It first requires to introduce the velocity 
fields [38] or the mass-current densities [39] in the fluid. The normal and superfluid 
components can then be extracted from the spatial correlation functions of this mass- 
current density [39]. Finally an analysis using renormalization group arguments leads 
to the so-called universal jump for the superfluid density: This density takes the value 
n s = 4/ A 2 on the low temperature side of the transition point, and n s = on the high 
temperature side. The existence and the value of the universal jump in the superfluid 
density was formally derived by Nelson and Kosterlitz [38] , and it was first confirmed to 
high accuracy in experiments with liquid He films by Bishop and Reppy [2]. 

This full description is outside the scope of this set of lectures. Here we simply 
illustrate how vortices drive the BKT transition, starting from a superfluid with finite 
n s , and then considering the free energy associated with spontaneous creation of a single 
free vortex. Without any loss of generality, in order to simplify the calculations we 
consider a circular geometry, with R — > oo the radius of the system. The kinetic energy 
cost of a single vortex placed at the origin is simply given by: 




where as always we assume that only the superfluid component rotates under the influence 
of the vortex. The normal component does not have any phase stiffness and its motion 
is not affected by the presence of the vortex. 

The entropy associated with a single vortex core is given by the number of distinct 
positions where a vortex of radius £ can be placed in a disc of radius R: 

(50) S = k B ln(^pj =2fc B ln 

Note that in the above calculations we ignore the "edge effects" such as the correction 




18 



Zoran Hadzibabic and Jean Dalibard 



to the energy for an off-centered vortex. One can check that these effects are negligible 
for R » £. Combining eqs. (49) and (50), we get for the free energy F = E - TS: 



We thus see that the free energy associated with a free vortex changes sign at n s A 2 = 4. 
Since ln(i?/£) diverges with the size of the system, this point separates two qualitatively 
different regimes. For n s X 2 > 4, F is very large and positive, so the supcrfluid is stable 
against spontaneous creation of a free vortex. On the other hand, for n s X 2 < 4, the large 
and negative F signals the instability against proliferation of free vortices. Appearance 
of first free vortices reduces n s and makes the appearance of further free vortices even 
easier, and this avalanche effect renormalizes the superfluid density to zero. We thus find 
that in contrast to 3d, where below the BEC critical temperature the superfluid density 
grows smoothly, in 2d the superfluid density cannot have any value between 4/A 2 and 0. 

Even though it does not explicitly address the microscopic origin of the vortex (the 
breaking of a vortex pair), this simple calculation correctly predicts the result for the 
universal jump in the superfluid density which takes place at the transition temperature 



This success relies on the fact that the above derivation is a powerful self-consistency 
argument. Whatever the origin of the vortex, and the relation between n s and total 
density n, it shows that it is inconsistent to suppose that we have a system with superfluid 
density which is non-zero, but smaller than 4/A 2 . Remarkably, this result also does not 
depend on the strength of interactions g, even though we know that the phase transition 
is mediated by interactions, since it does not occur in the ideal gas. Recent classical field 
Monte-Carlo calculations performed with parameters relevant for atomic gases [40, 41, 42] 
have confirmed the proliferation of vortices around the critical point characterized by 
eq. (52), although the transition was rounded off by finite-size effects. 

If we repeat the above arguments for tightly bound vortex pairs, we find that a finite 
density of pairs is present in the gas at any non-zero temperature. The energy of a pair is 
finite since the velocity field decays asucc 1/r 2 and the integral J v 2 d 2 r is convergent. 
On the other hand, the entropy is still divergent and essentially identical to the result of 
eq. (50), since the size of a tightly bound pair is of the same order as the size of a single 
vortex. The free energy for vortex pairs is therefore always negative. At any non-zero T 
pairs are continuously created and annihilated through thermal fluctuations. 

As the temperature is increased, but still kept below Tbkt, the density of pairs grows 
and also thermal fluctuations result in pairs of increasing size. As the average size of the 
pairs becomes comparable to the distance between the pairs, they start to overlap and 
this leads to effective screening of the attraction between two bound vortices, making 
it easier for fluctuating pairs to grow to even larger sizes. We can use an analogy with 



(51) 




(52) 



n s A 2 = 4 . 
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a Coulomb gas: Two nominally paired but well separated vortices create a field which 
polarizes the more tightly bound vortex dipoles between them. This results in an effective 
dielectric constant which reduces the attraction between two oppositely charged vortices. 
As Tbkt is approached from below, this creates an avalanche effect which eventually leads 
to breaking up of pairs and creation of a plasma of free vortices. Within the Coulomb 
gas analogy, this plasma provides a perfect screening at the transition point: A charge 
added at a given point does not change the flux of the electric field across a large radius 
circle centred on this point. Proper BKT calculation identifies the phase transition with 
the temperature at which the average size of the pairs, or equivalently the screening 
dielectric constant, diverges. 

3'3. Results of the microscopic theory. - The relation (52) between the superfluid 
density and the temperature Tbkt at the critical point is elegant and universal, in the 
sense that it does not depend on the interaction strength g. However this self-consistent 
result alone does not allow us to predict the value of Tbkt in a given system. It just 
tells us that whatever Tbkt is, the superfluid density jumps accordingly to 4/A 2 at the 
transition. 

Calculating the actual value of Tbkt in terms of the bare system properties n and g 
is a difficult problem, because it depends on the short-distance physics, such as density 
fluctuations which control the relationship between n s and n at the transition. In the 
weak coupling limit, j<l,a combination of analytical [43] and numerical [29, 44] efforts 
gives the value for the critical phase-space density: 

(53) D c = (n\ 2 ) c = \n(C/g) , 

where the dimensionless constant C = 380 ±3 is obtained by a classical field Monte-Carlo 
simulation [29]. The calculation leading to eq. (53) is formally valid only in the weak 
coupling limit g -C 1. We can set one obvious bound on its validity by noting that at 
the transition n > n s . Setting D c > 4 we obtain g < 7. This result is remarkably close 
to our estimate of the strong coupling limit g = 2n (eq. (22)). 

The numerical simulation of [29] provides another quantity of interest, which char- 
acterizes the reduction of density fluctuations due to interactions. The authors of [29] 
introduce the quasi- condensate density defined as: 

(54) n qc ^(2n 2 -(n 2 (r))) 1/2 . 

When interactions are negligible, (n 2 ) = 2(n) 2 and n qc = 0. On the other hand if density 
fluctuations are completely suppressed, (n 2 ) = (n) 2 and n qc = n. According to [29], at 
the critical point 

n qc 7.16 



(55) 



n HC/~g) 



This result indicates that n qc is of the order of the total density n at the transition 
point, unless the interaction strength g is exponentially small. In other words, for realistic 
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parameters density fluctuations are notably reduced in the vicinity of the BKT transition, 
which justifies the simplified Hamiltonian (26) used above. Actually this result sets a 
stronger constraint than eq. (53) on the applicability of the classical Monte-Carlo analysis: 
the condition n > n qc requires g < 0.3. 

Note that the terminology quasi- condensate density can sometimes be misleading. 
The quantity n qc takes a non-zero value even above the critical temperature for the 
BKT transition. It refers only to the properties of the density distribution in the gas, 
and not to the phase distribution as the word condensate might suggest. For example, 
we can have n qc <~ n at T > Tbkt, but this does not imply that a large contrast 
interference would be observed if one would superpose a pair of 2d gases prepared in the 
non-superfluid regime (see Section 7). 

Finally, we note that as the transition temperature is approached from above, the 
length scale I characterizing the exponential decay of correlations gi(r) <~ e~ r l it in the 
normal state diverges as: 



where a is a model-dependent dimcnsionless constant. Diverging correlation length is a 
very general property of phase transitions, but while in the case of most conventional (3d) 
second-order phase transitions the divergence of the correlation length is polynomial, in 
the case of the BKT transition it is exponential. This makes the critical region above 
Tbkt larger, and has implications for the broadening of the transition in finite size 
systems (see e.g. § 4'3). 

4. The 2d Bose gas in a finite box 

It is well known that finite size effects can play a significant role in the quantitative 
analysis of the phase transitions that are observed experimentally. In the two-dimensional 
situation of interest here, this is even more the case because the thermodynamic limit is 
reached only when ln(i?/£) >• 1 (see for example eq. (51)), which can be only marginally 
true for realistic systems. The analysis of these finite-size effects is therefore crucial for 
understanding the observed phenomena, and we will do so for the cases of a flat box 
potential (this section) and a harmonic trap (next section). 

We consider in this section a gas of N particles confined in a flat box of area L 2 . The 
box is supposed to be square (L x = L y = L) except in the last subsection (§ 4'5) where 
we discuss possible effects due to an anistropic confinement (L x ^ L y ). The confinement 
introduces a natural energy scale E = H 2 /(mL 2 ) in the problem and makes it possible 
to reach a true Bose-Einstein condensate at non-zero temperature, in contrast to the 
infinite case. In this section we first review the results that can be derived for the ideal 
gas case, and then discuss what happens for an interacting system. We will assume that 
the size L is much larger than the thermal wavelength A so that E <C k^T. 



(56) 
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4T. The ideal Bose gas. - In absence of interactions, the statistical description of a 
Bose gas in a square box of size L is straightforward. Let us choose for simplicity periodic 
boundary conditions so that the single particle eigenstates are plane waves e lk r /L, of 
energy tk — h 2 k 2 /2m, where the momentum k = (j x ,j y )(2ir/L), with j XtV positive, zero 
or negative integers. The chemical potential /i is always negative so that the fugacity 
Z = e 13 ^ lies in the interval < Z < 1. Three regimes can be identified, the first two 
being identical to what we have met for the infinite case: 

• The non-degenerate, high temperature regime with a phase space density D = 
nX 2 <C 1. This corresponds to a negative chemical potential such that \/j,\ ^> k^T 
(Z <C 1). The one-body correlation function is a gaussian function in this regime, 
gi(r) = ne~ 27Tr / x , and is vanishingly small for distances r ^> A, which are much 
smaller than the box size L. The confinement has no significant consequence on 
the coherence of the gas in this regime. 

• The degenerate, but non-condensed regime, where the momentum distribution is 
bimodal, with a Lorentzian shape for small k (kX -C \/4tt) and a gaussian shape for 
large k. The one-body correlation function decays exponentially at large r, with 
a characteristic decay length I = XeP I 2 j \fH . No significant condensed fraction 
appears as long as I is small compared to the size L of the sample, i.e. when 
D < ln(47ri 2 /A 2 ) or cquivalently |/x| > E /2. 

• The condensed regime, which occurs when the characteristic decay length I of g\ is 
larger than the system size L. This occurs when the phase space density D reaches 
the value ln(47rL 2 /A 2 ) (or equivalently \fi\ < E /2). A significant phase coherence 
then exists between any two points in the gas. 

4'2. The interacting case. - We now turn to the interacting case with repulsive in- 
teractions and discuss what can be expected in the vicinity of the BKT transition. For 
now we assume that the size of the sample is large enough so that D <C ln(47rL 2 /A 2 ) 
at the point where phase space density D is equal to D c , the critical phase space den- 
sity for the BKT transition in an infinite system (see eq. (53)). This condition has the 
following physical meaning: suppose that we increase the density of particles n at fixed 
temperature; the point where the superfluid transition occurs is reached well before the 
point at which the Bose-Einstein condensation due to the finite system size would occur 
in absence of interactions (see section 4'1 above). 

We first recall the nature of the superfluid transition in an infinite, homogenous 
sample. When D is notably below D c , but larger than 1, one expects that no superfluid 
component is present and <?i (r) decays exponentially. When D reaches D c the superfluid 
transition occurs and gi(r) decays algebraically (eq. (48)): gi(r) ~ n s (^/r) a for r > 
£, with a = l/(n s X 2 ) and n s A 2 > 4. According to the Penrose-Onsager criterion, 
no condensate is expected in an infinite system since gi(r) vanishes at infinity for any 
non-zero temperature. In sharp contrast with the infinite case, we now show that the 
BKT transition in a realistic finite system is always accompanied by the appearance of 
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a significant condensed fraction, denned as the largest eigenvalue n of the one-body 
density matrix. The basic reason for this effect is that the algebraic decay of gi(r) is 
extremely slow. To prove this result we proceed in two steps: first we give a general 
relation between IIo and the value of gi(r) for distances r comparable to the size L of 
the box; then we discuss what a realistic value of IIo can be for a typical atomic gas. 

Let us denote by Uj and <f>j (r) the eigenvalues and eigenstates of the one-body density 
matrix. The condensed fraction n is associated with the eigenstate 4>o( r ) = 1/L. We 
now consider the general expansion of g\ 

(57) Si(r) = AT^n^(0)^(r) 

3 

and integrate this expression over the area L 2 centered on the origin. For simplicity we 
integrate the left hand side over a disk of radius R = L/^/tt and the right hand side over 
a square of side L. This simplification cannot significantly affect our conclusions. The 
left hand side gives 

(58) J gi (r)d 2 r = 2irn s j\i/r) a r dr ~ L 2 ^— 9l (L) . 

On the right hand side, only the contribution of j = is non zero and gives ATIo. All 
the </>j's with j ^ are orthogonal to (f>o, so their integral over L 2 is zero. We thus get 
gi(L)/n <~ n . Using eq. (18) we can also write this result for the condensed fraction as 
IT ~ (n s /n)g~ a / 2 N- a / 2 . 

For a < 1/4, in practice we have g~ a ^ 2 <~ 1 and n s ~ n, so just below the transition 
temperature IIo ~ N^ 1 / 8 . Taking N = 10 5 as a typical value for cold atom experiments, 
we get n ~ 0.25. ( 6 ). We therefore meet here a paradoxical situation: the appearance 
of a non-zero condensed fraction may be used as a signature of the BKT transition, 
whereas the BKT mechanism was presented (for an infinite system) as a feature that 
takes place in a 2d interacting gas instead of the usual BEC of 3d Bosc fluids. Note that 
in a "true" BEC the condensed fraction n should not explicitly depend on N. However, 
for a < 1/4 this distinction becomes experimentally irrelevant, and in order to observe a 
BKT transition with no significant BEC one would need to consider unrealistically large 
samples. This was pointed out by the authors of [45] who wrote the famous statement (in 
the context of 2d magnetism): 'With a magnetization at the BKT critical point smaller 
than 0.01 as a reasonable estimate for the thermodynamic limit, the sample would need 
to be bigger than the state of Texas for the Mermin-Wagner theorem to be relevant!". 
A similar remark holds in the context of superfluid helium films [2] . 



( 6 ) We get an equivalent estimate from YI ~ g 1 (L)/n and n s ~ n in the algebraic decay regime. 
In cold atom gases, the typical values of A and £ arc 0.1 — 1 fim, while the maximal system size 
is L ~ 100 urn. At the transition point we have II ~ (10~ 3 ) 1/4 to (10~ 2 ) 1/4 ~ 0.2 to 0.3. 
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4'3. Width of the critical region and crossover. - The intricate mixing between the 
BKT mechanism and the emergence of a significant degree of coherence exists even at 
temperatures slightly above the BKT transition point. We mentioned in § 3'3 that in 
the normal state the characteristic decay length £ of gi(r) diverges exponentially in the 
vicinity of the critical point: \n(£/X) w (o,Tbkt/(T — Tbkt)) 1//2 , where a is a model- 
dependent coefficient (see eq. (56)). The critical region where I becomes larger than A as 
a precursor of the BKT transition is therefore very broad, [T — Tbkt) ~ Tbkt- Further, 
there clearly exists a temperature close to (but still above) Tbkt for which £ exceeds the 
system size. At this temperature a significant condensed fraction appears in the system. 
Because of the exponential variation of I with T — Tbkt, the temperature range where 
I > L, and condensation gradually sets in, can be significant ( 7 ): 

AT _ AD a 
Tbkt" ~ ^7 ~ (m(T/A)) 2 " 

Although cold atom systems are so far commonly confined in harmonic rather than 
box-like potentials, it is interesting to provide an estimate for this case. Taking a = 1 
and reasonable values for L/X between 10 and 100, the BKT transition is expected to 
become a crossover with a relative width AT/T B kt ranging from 5% to 20%. 

The above analysis explicitly concerns only the emergence of a non-zero condensed 
fraction, but it also suggests broadening of the universal jump in the supcrfluid density. 
For example, using standard Bogoliubov argument (see § 2'4) we can deduce that finite 
condensed fraction in an interacting system also implies a finite superfluid density. Also, 
if we take the definition of superfluid density which associates it with the energy cost 
of twisting the phase of the wave function at the edge of the system [46, 47], we again 
conclude that £ > L implies a finite superfluid density. In the critical region, quantitative 
conclusions might actually depend on what theoretical definition of superfluid density we 
accept, but the qualitative conclusions will not change. 

4'4. What comes first: BEC or BKT? . - This is an often raised and subtle question. 
We have so far discussed the case of a large system such that D c <C ln(47rL 2 /A 2 ) where 
D c = ln(380/<?) is the critical phase space density for the BKT transition in an infinite 
system. For such large systems, the first relevant mechanism that occurs when increasing 
the phase space density is a BKT transition. However we have seen that the approach of 
the BKT threshold always results in the appearance of a significant condensed fraction. 
We can qualify this "BKT-driven" condensation as "interaction-enhanced" , since it would 
not take place in an ideal gas with the same density and temperature. Experimentally, 
the strength of interactions in cold atom systems can be dynamically controlled using a 
Fcshbach resonance [48, 49]. One can therefore imagine preparing a non-interacting gas 
at phase space density where no condensation occurs and then driving the condensation 



( 7 ) Here we can assume that A under the logarithm on the right hand side is constant over the 
range AT. 



(59) 
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via the BKT mechanism by turning on the interactions. In the opposite regime of a small 
system where the BKT transition would require D c > ln(47rL 2 /A 2 ), the first phenomenon 
that is encountered when the phase space density is increased is "conventional" Bose- 
Einstein condensation as for an ideal Bose gas. As in the 3d case, in presence of weak 
repulsive interactions the formation of a condensate is accompanied by the apparition of 
a superfluid fraction with a comparable value. 

It would be very interesting to study cold atomic gases in a (quasi-)uniform potential 
and vary the experimental parameters so as to explore both regimes. However we point 
out that this will not be an easy task if we require that the criteria for the two types of 
transitions are well separated, for example by more than the crossover width discussed 
in § 4'3. To analyze the system requirements for reaching the two different regimes, it is 
convenient to fix the ratio L/X and write ln(47rL 2 /A 2 ) = -fD c . Now 7 is a dimensionless 
parameter such that 7 > 1 means that (as particle number is increased) condensation 
occurs via the BKT mechanism. The critical number for the BKT transition is then 

(60) N C = ^D C = ^D c e> D < , 

while the critical number for condensation in an ideal gas is "fN c . For illustration purposes 
we may define the BKT regime by 7 > 1.5, and the BEC regime by 7 < 0.5. (For values of 
7 close to 1, the two effects are difficult to disentangle experimentally.) The experiments 
with cold atoms have so far been performed at coupling strength g ~ 10~ 2 — 10 _1 . Taking 
g = 0.1, 7 = 1 corresponds to iV c ss 2.5 x 10 3 , and the BKT regime 7 = 1.5 corresponds 
to iV c w 1.5 x 10 5 , which is easily achievable. However the opposite BEC regime of 
7 = 0.5 corresponds to N c w 40; studying such a small particle number is experimentally 
very challenging, although it might become feasible with the development of single-atom 
detection [50, 51, 52]. For a more weakly interacting gas with g = 0.01, 7 = 1 corresponds 
to N c w 3 x 10 4 , and 7 = 0.5 to iV c w 160, which might be easier to explore. On the other 
hand 7 = 1.5 corresponds to iV c w 6 x 10 6 , which would be experimentally challenging. 
It would therefore generally be difficult to explore both the large (BKT) and the small 
(BEC) system regime using the same value of g, and reaching the BEC regime may 
require a more weakly interacting quasi-2d atomic gas than has so far been studied. 

4'5. The case of anisotropic samples. - So far we have assumed that for phase space 
densities D larger than the critical value D c for the BKT transition, the functional form 
of g\ found in the infinite case (algebraic decay) remained valid for a finite size system. 
This assumption is reasonable for square samples (L x — L y ), but may not be valid for 
anisotropic samples, with a width along one direction (say x) much larger than the other 
one: L x ^> L y . We briefly review the expected properties in this regime, which is relevant 
for several of the previous or current experimental setups. Since we are interested in the 
regime D > D Cl we assume that a superfluid component is present in the sample and 
we use the hamiltonian Hg given in (41) to estimate the amplitude of phase fluctuations 
and their consequence on the one-body correlation function g\. 
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We start from the result (46) obtained in the infinite case. In a finite size system 
the integral over k is replaced by a discrete sum over k — 2ir(j x /L x , j y / L y ) times the 
constant prefactor 4ir 2 / (L x L y ): 

, . , . , ., , 2-7T v-^ 1 — cos(fc • r) 

(61) ^ (r)/ ^,____£_J_J. 

We are interested here in the decay of g\ along the long axis of the sample and we choose 
r = xu x , where u x is the unit vector along the x direction. We take x <C L x so that 
the finitencss of the sample along x has no relevance here. On the contrary the size L y 
(<C L x ) will play an important role since we can choose x either small or large compared 
to L y . We now show that these two cases lead to different decaying regimes for g\. 

Since we assume \x\ -C L x , the discrete sum over k x can always be replaced by an 
integral and we get: 

(62) InMaO/n.) - ' v , 1 -«»(*,*) , t _ * ^1-e"^ 



1 \ - f 1 - cos{k x x) tt y" 

~n s \ 2 L y ^J kl + k 2 y dkx ~ n s \ 2 L y ^ \k y \ 



We single out the contribution of k y = in the sum, introduce a cutoff fc max for large k y 
y™i Vi » Injmax and Y.% 



and use £'=T 1/j w In j max and CVi = - ln(l - C) for < C < 1- We then get 



(63) 



H^)/n s ) = + in (l - e- N^)] } 



Two regimes clearly appear in this expression. If ir\x\ <C L y then the logarithm on the 
right hand side is the dominant term, and we recover the algebraic decay that holds for 
an infinite system: 

n 1 
(64) fc"^ « < Ly : gl (x) » s with a 



Jxlfcmax)" Tl s X 2 

This result is intuitive: as long as we probe the coherence of the system on a distance 
shorter than the smallest size of the sample, the anisotropy introduces no significant 
deviation with respect to an infinite system. 

The situation is dramatically different when \x\ » L y . In this case the dominant 
contribution on the right hand side of (63) is the linear term n\x\/L y . This term, which 
originates from the contribution of the k y = mode to the sum (61), leads to an expo- 
nential decay of g\ : 

(65) L y < |a;| : gi(x) w n s e ^ xWd with d = n s \ 2 L v /ir . 

It is quite remarkable that when we probe the coherence of this anisotropic system on 
distances x > d > L y , we obtain an exponential decay as if the system was not superfluid. 
At the same time, the characteristic distance over which g\ decays, d, explicitly depends 
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on the superfluid density. The physical interpretation of this counterintuitive result is 
that over such distances the system acquires a quasi-one dimensional character: the phase 
stiffness between the origin and the point at coordinate x is decreased with respect to 
an infinite plane because there is a severe reduction in the number of independent paths 
connecting these two points. Ultimately, for very large |x|, only the channel k y = 
contributes significantly to the connection between these two points. This explains why, 
although the system is superfluid, the decay of g\ turns to an exponentially decaying 
function, that is characteristic of Id degenerate gases (see also [53] for a similar discussion 
for elongated atomic 3d gases). 

5. The 2d Bose gas in a harmonic trap 

Up to now, experiments performed with (quasi-)2d atomic gases used harmonic trap- 
ping in the xy plane. We discuss in this section how the presence of the harmonic trapping 
potential modifies the above results. We shall see that it can lead to a dramatic change of 
the properties of the system, through the modification of the density of states of the single 
particle Hamiltonian. In particular, in this case "conventional" Bose-Einstein condensa- 
tion, as defined through the saturation of excited states at some non-zero temperature 
T c , can occur in the ideal Bose gas even in the thermodynamic limit. However, in pres- 
ence of repulsive interactions, and in large enough systems, this type of condensation is 
suppressed and replaced by the BKT normal to superfluid transition. 

5T. The ideal case. - Consider for simplicity an isotropic 2d harmonic potential 
V(r) = muj 2 r 2 /2. The single particle energy levels are Ej = (j + l)fuv, with j positive 
or zero integer, and each level having a degeneracy gj = j + 1. The maximum number 
of atoms N c that can be placed in all excited states (j > 0) at a given temperature T is 
obtained by choosing a chemical potential \x equal to the ground state energy: 



where £ = huj/(kBT). Assuming £ -C 1, the discrete sum can be replaced by an integral 
and one obtains 



For an atom number N > N C (T) there must be at least N — N c atoms occupying the 
single particle ground state j — 0. Equivalcntly, for a given atom number TV placed in 
the trap, there must be a significant fraction of the atoms that occupy the ground state 
j = if the temperature is reduced below the critical value 



(66) 




(67) 




(68) 



k B T c = —hioVN 

7T 
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Since the ground state is separated from the first excited state by a non-zero gap hu>, 
this Bose-Einstcin condensation can be viewed as a natural consequence of the finite size 
of the system [54], similar to the condensation of the ideal gas in a finite box. However, 
one can see that the condensation of the ideal gas in a 2d harmonic trap is a more 
interesting phenomenon by considering the appropriately defined thermodynamic limit 
for the harmonic confinement. This limit is obtained by taking N — > oo and uj — > 0, 
while keeping T and Nui 2 constant. For a gas described by Boltzmann statistics, this 
ensures that the central density n = N mu 2 / {2tt k^T) remains constant. Eq. (68) 
leads to a non-zero critical temperature in the thermodynamic limit, contrarily to what 
happens in the uniform case. This result can be understood by noticing that the density 
of states has a different functional form in a box (p(E) constant) and in a 2d harmonic 
potential (p(E) oc E) [55]. The vanishing density of states at E = for a 2d harmonic 
potential leads to a similar situation to the 3d uniform case, hence the possibility for a 
genuine Bose-Einstein condensation in the ideal gas (for a discussion of small logarithmic 
anomalies in the compressibility of the trapped ideal 2d gas, see [56]). 

Quite remarkably the result (67) can be recovered by starting from the uniform result 
given in cq. (6): n = — A~ 2 ln(l — e 13 ^) and using a local density approximation (LDA). 
This approximation amounts to replacing the uniform chemical potential p, by the local 
one n — V(r), which gives the following expression for the total atom number: 



where we set R = r/r? with rip = k^T/rruJ 2 . For n = 0, the result coincides with 
cq. (67) . Therefore in spite of the fact that LDA leads to a diverging spatial density at 
the center of the trap for p = [n(r) oc — ln(r)], it provides the same upper bound N c 
as eq. (66) for the total number of atoms, assuming no macroscopic occupation of the 
single particle ground state. 

5'2. LDA for an interacting gas. - In order to take into account repulsive interactions 
for a trapped gas, we use again the local density approximation. We will start with an 
analytical mean-field treatment based on the Hartree-Fock approximation. We will then 
use the numerical results of a classical field Monte-Carlo approach [44] that will provide 
a more precise determination of the BKT transition. 

In the mean- field Hartree-Fock approach when no condensate is present, interactions 
are taken into account by adding the energy 2gn(r) to the external potential [57, 58]. 
The local chemical potential is now p — V(r) — 2gn(r) so that the local phase space 
density D(r) = n(r)X 2 is the solution of the implicit equation 



N =-X 



-2 




(69) 




(70) 



D(r) = -ln{l - Zcxp[-(3V{r) - gD(r)/n]} 
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Putting R = r/rx as above, we can write the total atom number as 



(71) 



N 



(id) 



6 f + °° 

-2 / D(R) RdR 



Jo 



where D is solution of 



(72) 



D(R) 



ln{l 



Zexp [-R 2 /2-gD(R)/n 



•]} • 



Note that the solution D(R) depends only on the fugacity Z and the interaction strength 
g. The trap frequency and the temperature do not appear explicitly so that the scaling 
of the atom number with co and T (at fixed Z) is identical to the result (69) for the ideal 
gas. 

Interactions, when treated at the mean-field level, dramatically change the nature 
of the solution of cqs. (71)-(72). For a given trapping frequency to and temperature 
T, and for any non-zero g, the atom number N obtained from eq. (71) can be made 
arbitrarily large by choosing properly the fugacity Z. The condensation phenomenon 
that was obtained in the ideal gas case does not occur anymore. This can be understood 
qualitatively. For an ideal gas, the saturation of the atom number occurs when the 
central density in the trap becomes infinite. In presence of repulsive interactions, this 
singular point cannot be reached and the mean-field treatment provides a solution for 
any atom number [59]. 

We have plotted in figure 2a the prediction of the mean-field approach for the central 
phase space density D(0) as a function of the total number of atoms in the trap, for 
various values of g. As expected, D(0) is a monotonically increasing function of the 
atom number TV: more atoms in the trap lead to a larger central density. The other 
expected feature is that, for a given atom number, D(0) decreases as the repulsion 
between atoms is increased. In particular the divergence of D(0) that is found in the ideal 
case for N = Nc^ does not show up anymore in presence of repulsive interactions. One 
could try to push further the mean-field analysis of the equilibrium state, and look for 
dynamical or thermodynamical instabilities that could appear above some critical atom 
number [60, 61, 62, 63]. However we will rather follow the spirit of LDA and assume 
that the normal to supcrfluid BKT transition occurs at the center of the trap when the 
phase space density at this point exceeds the critical value D c (eq. (53)) predicted for 
the uniform system [64]. 

Within the mean-field Hartree-Fock analysis, one can show that, to a very good 
approximation, the number of atoms that have to be placed in the trap so that the 
central phase space density reaches the critical value D c is [65] 



(73) 




An obvious consequence of this result is that for a given trap and a given temperature, 
the BKT threshold in presence of interactions requires a larger atom number than the 
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Fig. 2. - Variation of the central phase space density as a function of the atom number, normal- 
ized by the critical atom number for the ideal gas in the same potential and at the same 
T. The value of the interaction strength g is given for each curve. The black squares indicate 
the values of N at which the BKT criterion of eq. (53) is met at the center of the trap, (a) 
Results obtained using the mean-field Hartree-Fock approach, (b) Results obtained using the 
bulk results of [44] and the local density approximation. 



BEC of the ideal gas. Equivalently, for a given atom number, the superfluid transition 
temperature in presence of interaction is lower than the ideal gas condensation temper- 
ature. In figure 2a we have indicated with black squares the value of Nc" 1 ^ for various 
interaction strengths. 

So far we have relied on the mean-field approximation to obtain the relationship 
between the density n(r) and the local chemical potential fx — V(r) — 2gn(r). Although 
this gives a good feeling for the scaling laws that appear in the problem, it cannot 
provide a very accurate description of the transition. Indeed the mean-field expression 
2gn(r) for the interaction energy can only be valid at relatively low density, where the 
density fluctuations are important, so that (n 2 ) = 2(n) 2 . When the density increases 
and/or the temperature decreases, density fluctuations are reduced and one eventually 
reaches a situation at very low temperature where those fluctuations are frozen out and 
(n 2 ) = (n) 2 . At zero temperature, one expects a quasi-pure condensate in the trap 
with a density profile given by the Thomas-Fermi law gn{r) = /i — V(r) (whereas the 
Hartree-Fock approximation would lead to replacing g by 2g in this equation). 

To capture the reduction of density fluctuations as the phase space density increases, 
we now use the numerical results of [44] obtained using a classical field Monte-Carlo anal- 
ysis. They provide the value of the phase space density D as a function of [ijkT in the 
vicinity of the BKT critical point for a uniform system. Injecting this numerical predic- 
tion in the LDA scheme, we obtain the results shown in fig. 2b for the central density as a 
function of the total atom number. As expected this figure is qualitatively similar to the 
one obtained using the mean-field Hartree-Fock approach. However the classical Monte- 
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Carlo results lead to a noticeable reduction of the critical atom number with respect to 
the mean field treatment. For example, for g = 0.15 (as in the ENS experiment [27], 
see below) the critical atom number for reaching the BKT threshold is expected to be 
~ 1.47V c (id) using the numerical predictions of [44] instead of ivi™^ ~ 1.9 A^ ld ^ using the 
Hartree-Fock approximation. One might wonder if the classical Monte-Carlo simulations 
of [44] , which assume g <C 1 , remain accurate for the relatively large interaction strength 
g = 0.15. The (positive) answer was given in [96], which provides a detailed comparison 
between the predictions of [44] and those of a quantum Monte-Carlo simulation of an 
assembly of trapped bosons for the interaction strength and trapping geometry of the 
ENS experiment. 

5'3. What comes first: BEG or BKT? . - In the previous section devoted to the study 
of a square potential we have explained that the distinction between a conventional BEC 
transition and a BKT transition is subtle, and we introduced two related concepts: 

• "BKT-driven condensation", meaning that if the phase space density is increased 
at constant g the first many-body mechanism encountered is the BKT transition, 
but due to the resulting slow decay of gi this transition is accompanied by the 
appearance of a finite condensed fraction. 

• "Interaction-enhanced condensation", meaning that there exists a range of phase 
space densities for which no condensation occurs in an ideal gas but condensation 
via the BKT mechanism can be induced by increasing the interaction strength from 
to ~g. 

In the case of a box potential these two concepts are equivalent. We identified a range of 
parameters, such that D c <C ln(47rL 2 /A 2 ), for which both effects occur. In the opposite 
regime of a small system and/or small g neither of the two effects occurs. 

The question "what comes first" is even more subtle in the case of a harmonically 
trapped gas because of the inhomogeneous density profile. In this case the notions of 
BKT-driven and interaction-enhanced condensation are not equivalent, and the answer 
depends on what we keep constant in an experiment, i.e. which path we follow in the 
phase diagram. Since in a harmonic trap ideal gas condensation occurs even in the ther- 
modynamic limit (i.e. if we neglect the discreteness of single-particle energy levels in the 
trap), we start by analyzing that case, and separately consider two different experimental 
paths: 

• The critical phase space density D c for a BKT transition (at a fixed non-zero 
g) is finite, while the critical phase space density for ideal gas condensation is 
infinite. Therefore in a standard experiment where g is kept constant and phase 
space density is increased, BKT-driven condensation always occurs. In this sense, 
in practice "BKT always comes first" . 

• While the critical phase space density for the BKT transition (at non-zero g) is lower 
than for the ideal gas condensation, the critical atom number at fixed T is higher. 
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The critical atom numbers iV c and iV c scale similarly with the temperature 
and the trap frequency (oc (k B T/huj) 2 ), and the ratio ivi BKT '/A r i ld ^ is always larger 
than 1. This can be seen from the mean-field result of eq. (73) or from the Monte- 
Carlo data shown in figure 2b. This means that, at fixed N, interactions always 
reduce the transition temperature. Therefore, contrary to the case of a square box 
potential, we can never have interaction-enhanced condensation. 

Note that it is not inconsistent that the BKT transition occurs at a lower critical 
density but higher critical number than the ideal gas BEC, because in a harmonic trap 
with fixed N and T the peak (phase space) density in a repulsively interacting gas is lower 
than in an ideal gas. Also note that if we work within the BKT theory and then formally 
take the g — > limit, we exactly recover the criterion for ideal gas condensation, which 
is usually derived from a conceptually completely different viewpoint of the saturation of 
single-particle excited states. We can therefore think of the BEC transition as a special 
non-interacting limit of the more general BKT theory. This connection naturally emerges 
when analyzing the case of a harmonically trapped gas, but it could not be made in a 
uniform system, where the critical temperature for both transitions vanishes in the 54O 
limit. 

Finally, we briefly comment on the case of a realistic experimental harmonic trap, 
where the spacing of the single-particle energy levels is non-zero. The results for the 
critical atom numbers for the BKT and the ideal gas BEC transition are essentially 
unaffected by the non-zero level-spacing. It therefore remains true that interaction- 
enhanced condensation is not possible. However, the ideal gas BEC in this case occurs 
at a finite phase space density Z?bec m the trap center, which can in principle be lower 
than D c for some values of g. In this case BKT-driven condensation would also no longer 
occur, and "BEC would come first" no matter what path we take in the phase diagram. 
This scenario is however not relevant for the currently realistic experiments. The value 
of -Dbec is not universal and depends on the details of the trapping potential, but we 
have evaluated it numerically for a typical trap used in the ENS experiments [26, 27], 
and obtained Dbec ~ 13 [66]. This means that the condition Dbec < D c can be fulfilled 
only in an extremely weakly interacting gas with g < I0~ 3 . Experimentally, this regime 
is essentially indistinguishable from the g — > limit, where the BKT and the BEC 
transition are no longer distinct. 

5'4. Width of the crossover. - The divergence of the correlation length that we already 
discussed in the case of a square box potential (see § 4'3), must also be taken into account. 
Suppose that one lowers the temperature of a cold gas in a trap until the BKT threshold 
density is reached right at the center of the trap. Since the density is everywhere lower 
than the critical threshold for BKT, LDA would imply that no significant supcrfluid 
fraction is present in the gas at this stage. However a significant part of the gas may 
exhibit a certain degree of coherence as we show now. The critical length (56) is now 
position dependent, since the critical temperature is a function of the local density. 
Since l(r) is a monotonically decaying function of the distance from the trap center r, 
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we can self-consistently assume that the gas is coherent over a region of radius r c such 
that £(r c ) — r c . We can also provide a crude estimate of r c : r c w r T (ln(r T /A)) -1 / 2 . 
For practical parameters (r-r/A ~ 10 — 100), we find r c ~ r-r, which means that this 
coherence actually extends over a significant fraction of the cloud when D = D c at the 
center of the trap. We can also estimate the width of the cross-over over which the 
condensed fraction becomes significant. If instead of taking D = D c at the center of 
the trap, we take D = 0.7 D c then I w 5A at the center of the cloud and no significant 
coherence exists at this point. The above analysis is confirmed at least qualitatively 
by numerical simulations performed using a classical field Monte-Carlo analysis. These 
simulations indeed indicate the emergence of an extended coherence over the cloud at 
temperatures 10% to 20% above the one for which the bulk BKT criterion is met at the 
trap center [42]. 

6. Achieving a quasi-2d gas with cold atoms 

The experimental realization of a 2d atomic Bose gas is based on a strongly anisotropic 
trap with one very tightly confining direction, say z, and two more loosely confined 
degrees of freedom, x and y. The z degree of freedom can be considered as frozen from 
the thermodynamic point of view if the energy gap A z between the ground state and 
the first excited state of the z motion is much larger than both k^T and the interaction 
energy gn (both being typically on the order of one to a few kHz) . Since the confinement 
along z is usually harmonic, with frequency co z , the gap is A z — Huj z . The z degree 
of freedom is thermodynamically frozen when the extension of the ground state of the 
z-motion, a z — yjh/muj z , is such that a z <C £, X/V2w. 

6T. Experimental implementations. - Conceptually, the simplest scheme to produce 
a 2d gas is to use a single gaussian light beam that is red-detuned with respect to the 
atomic resonance. The beam propagates along the x direction, with waists along y and 
z such that w y >• w z , so that it forms a horizontal light sheet. The dipole potential 
created by this light sheet attracts the atoms towards the focal point, and ensures a 
strong confinement in the z direction. This technique was used at MIT to produce the 
first atomic gas (of sodium atoms) in a quasi-2d regime [67]. More recently it has been 
implemented at NIST to study the coherence properties of the 2d gas [28] . 

One can also produce a 2d gas by using an evanescent light wave at the surface of a 
glass prism [68, 69], so that the atoms are trapped at a distance of a few micrometers from 
the horizontal glass surface. The confinement in the horizontal xy plane is provided by an 
additional laser beam or by a magnetic field gradient. The fact that the confinements in 
the xy plane and along the z axis have different origins is an interesting feature because 
it offers the possibility, by releasing only the planar confinement, to study the ballistic 
expansion of the atoms in the xy plane only (see § 7'2). Another experimental system 
providing independent confinement in the xy plane and along z has been investigated at 
Oxford, where a blue detuned, single node, Hermite Gaussian laser beam traps atoms 
along the z direction, and the confinement in the xy plane is provided by a magnetic 
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field gradient [70]. 

Two-dimensional confining potentials that are not based on light beams have also 
been investigated. One possibility discussed in [71] consists in trapping paramagnetic 
atoms just above the surface of a magnetized material that produces an exponentially 
decaying field. The advantage of this technique lies in the very large achievable frequency 
uj z , typically in the MHz range. One drawback is that the optical access in the vicinity 
of the magnetic material is not as good as with optically generated trapping potentials. 
Another appealing technique to produce a single 2d sheet of atoms uses the so-called 
radio-frequency dressed state potentials [72, 73, 74]. 

A Id optical lattice setup, formed by the superposition of two running laser waves, is 
a very convenient way to prepare stacks of 2d gases [75, 76, 77, 78, 79]. The Id lattice 
provides a periodic potential along z with an oscillation frequency lo z that can easily 
exceed the typical scale for interaction energy and temperature. The simplest lattice 
geometry is formed by two counter-propagating laser waves, and provides the largest lo z 
for a given laser intensity. One drawback is that the lattice period is small (Ao/2, where Ao 
is the laser wavelength) so that many planes are generally populated and the addressing 
of a single plane is difficult. Therefore practical measurements only provide averaged 
quantities. Another interesting geometry consists in forming a lattice with two beams 
crossing at an angle 9 smaller than 180° [80]. In this case the distance A /(2 sin(0/2)) 
between adjacent planes is adjustable, and each plane can be individually addressable if 
this distance is large enough [50, 81]. Furthermore the tunneling matrix element between 
planes can be made completely negligible, which is important if one wants to achieve a 
truly 2d geometry and not a periodically modulated 3d system. 

Finally, while here we are primarily interested in continuous 2d gases of spinless 
Bosons, two related experiments on 2d physics also need to be mentioned: 

First, an experiment performed in Boulder constitutes a direct implementation of 
the XY model [82]. There, an array of parallel elongated (quasi-)condensates is created 
in a 2d optical lattice, and tunneling matrix element J provides a Josephson-type cou- 
pling between the neighboring lattice sites. In this system proliferation of vortices is 
observed when the temperature is increased. Vortices are detected by turning off the 
optical lattice and allowing the quasi-condensates trapped on different sites to overlap 
and interfere. The measured surface density of vortices as a function of the ratio J/T is 
in good agreement with the BKT theory applied to this system. 

Second, in an experiment at Berkeley 2d physics was studied in a spinor Bose-Einstein 
condensate of Rb atoms with total spin F — 1 and weak ferromagnetic spin-dependent 
interactions [83]. This system is anisotropic, but still 3d with respect to the density 
degrees of freedom, i.e. the healing length £ is shorter than the shortest extension of 
the cloud, along z. However, weak spin-dependent interactions correspond to a longer 
healing length £ s , so that the system is 2d with respect to the spin degrees of freedom. 
In this case the magnetization transverse to the quantization axis has a role analogous to 
the phase of the wave function in a spinless Bose gas. At low T ferromagnetic interactions 
favor spontaneous symmetry breaking but spin-vortex structures are also observed. 
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6'2. Interactions in a 2d atomic gas. - To address the role of interactions in these 
gases, we start with some considerations concerning the quantum scattering of two atoms 
when the z motion is strongly confined. In a strictly 2d problem and at low energy, the 
scattering state between two identical bosonic particles with relative wave vector k is 
[25] 

/ j pikr a 

(74) Mr) ^ e ^r_ f{k) f{k) 



8n y/kr ' — ln^a 2 .) + in 

where 02 is the 2d scattering length. One should notice that contrarily to the 3d case, 
the scattering amplitude f(k) does not tend to a non-zero finite value when k tends to 
0. In the experimental implementations of 2d gases that have been achieved so far, the 
confinement along z was still relatively weak from a collisional point of view, in the sense 
that the thickness a z of the gas remained notably larger than the 3d scattering length 
a s . The scattering problem in this confined geometry has been discussed in [84, 85] (see 
also [86]); the general expression (74) for the scattering state V'fc remains valid and the 
scattering amplitude can be written 

(75) f(k) ( 1 



2na z /a s — ui(k k 2 a z ) + in 
with k ~ 3.5, corresponding to the 2d scattering length 



(76) a-i = a z \fn exp 



n a, 
2 a f 



Now for all experiments realized so far, the first term \f2na z /a s in the denominator of 
cq. (75) is large compared to 1, and dominates over the logarithmic term \n{K,k 2 a 2 z ) and 
the imaginary term in. We can then take a constant scattering amplitude (as in 3d) to 
describe the collisions in the gas: f(k) = g ~ \f&na s j a z . With this approximation the 
interaction energy of the gas with density n(r) in the xy plane is 

17, 

d'r . 



(77, ft. = |f / „>, 



The corresponding values for the 2d scattering length are extremely small, due to the 
exponential factor in cq. (76). Taking for example a z = 200 nm and a s = 5 nm ( 87 Rb 
atoms), we find = 6 10~ 29 m. Typical surface densities are in the range I0 13 m~ 2 , 
and the dimensionless parameter na\ that is relevant for perturbative expansions of the 
equation of state of the 2d Bose gas [87, 88, 89, 90, 91, 92, 93] is also extremely small: 
na\ ~ 4 10 -44 for the numbers given above. 

The expression (77) can also be obtained by starting from the 3d interaction energy 

(78) E intM = [ „§(r) al\ , 

771 J 
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in which we plug directly n 3 (x, y, z) — n(x, y) exp(— z 2 / a 2 z ) / ^Jita 2 . However the validity 
condition a s <C a z remains hidden in this procedure. 

To summarize, the collision dynamics in the experiments performed so far is still 
dominated by 3d physics. The 3d scattering length a s is much smaller than the thickness 
of the gas and the scattering amplitude is nearly fc-independent. This regime is often 
referred to as 'quasi-2d'. It is important to note that the term 'quasi-2d' is also used 
to describe another aspect of the 2d gases: very often the temperature of the gas (and 
possibly the interaction energy) is not small compared to Tvjj z , but comparable or even a 
bit larger. We discuss in the next subsection how to handle this problem. 

6'3. Residual excitation of the z-degree of freedom. - For a quantitative analysis of 
experiments performed with 2d gases, in particular for the determination of the tem- 
perature, it is important to take into account the residual excitation of the ^-degree of 
freedom. This was first pointed out in [94], where a quantum Monte-Carlo simulation 
gave an estimate for the distortion of the density profile due to this residual excita- 
tion. Several possible ways were subsequently proposed to take this excitation into ac- 
count [94, 65, 66, 95, 96]. The simplest method consists in renormalizing the interaction 
strength g to account for the density profile of the gas along the z direction [94]. The 
predictions derived with this method were compared with Quantum Monte-Carlo results 
in [65] and later analyzed in detail in [95]. In the following we outline the slightly more 
elaborate treatment of [66] which has the advantage of taking into account not only the 
thermal excitation of the z degree of freedom, but also the possible deformation of the 
ground state of the z-motion due to atomic interactions. 

The method used in [66] is a direct implementation of the Hartree-Fock approximation 
(see e.g. [97]) and we first present it for a gas which is uniform in the xy plane. We choose 
a 3d trial density profile n 3 (z) uniform in the xy plane and varying along the strongly 
confined z direction. We then consider the Hamiltonian with the mean-field energy 

(79) H = -^-V 2 + \rmj 2 z z 2 + 2g^n 3 (z) . 

The single particle cigenfunctions of this Hamiltonian can be written tp k ,j{x, y, z) = 
ipj{z) e^^+M) /2tt, with energy E kJ = h 2 k 2 /(2m) + £j , where k 2 = k 2 x + k 2 y . The 
normalized functions fj{z) and the energies ej of the z-motion of course depend on 
the choice of the trial density profile n 3 (z). In the Hartree-Fock approximation the 
average occupation of the single particle level ip k ,j is given by the Bose factor f(E k j) = 
(exp((3(Ei~,j — /«)) — We calculate the corresponding 3d density profile, which is 

still uniform in xy and has the following z dependence: 

(80) = £ I d 2 k\^\ 2 f(E kij ) = -^J2\^)\ 2 ^^- Ze~^) . 

j j 

The self-consistency of the Hartree-Fock approximation requires that n 3 (z) and n' 3 (z) 
coincide, which can be achieved by iterating the solution of the above set of equations 
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until a fixed point is reached. With this method, we fulfill two goals: (i) We take into 
account the residual thermal excitation of the levels in the z direction, (ii) Even at 
zero temperature we take into account the deformation of the z ground state due to 
interactions. When interactions can be neglected, the eigenstates ifij(z) are the Hermite 
functions and €j = hui z (j + 1/2). 

Using the local density approximation, the above method can be straightforwardly 
adapted to the case where a trapping potential V± is present in the xy plane. The trial 
density distribution n 3 (r) is now a function of all three spatial coordinates. At any point 
(x,y), we treat quantum mechanically the z motion and solve the eigenvalue problem for 
the z variable 



(81) 



-h 2 d 2 

2m dz 2 



+ V cS (r) 



(pj(z\x,y) = Cj(x,y) ip {z\x,y) 



where V e s(r) = V±(x, y) + muj 2 z z 2 /2 + 2g ( - 3d ^n 3 (r) and / \tpj(z\x, y)\ 2 dz = 1. Treating 
semi-classically the xy degrees of freedom, we obtain a new spatial density 

(82) ^(r) = -^^|^-(z|x,y)| 2 ln(l-Z e -^^)) . 

j 

Again the Hartree-Fock prediction is obtained by iterating this calculation until the 
spatial density n 3 (r) reaches a fixed point. Since ipj is a normalized function of z at any 
point {x,y), the total 2d density is 

(83) n(x, y) = J n 3 (r) dz = ~ (l - Ze~^ x ^) . 

j 

In the limit where only the ground state j = of the z motion is populated, the result of 
this Hartree-Fock approach coincides with the solution of (70). The method used in [95] 
is similar to this approach, but the deformation of the eigenstates due to mean-field 
interaction was neglected. 

The density profiles predicted by this method have been compared with the results 
of a quantum Monte-Carlo simulation [96]. An important outcome for the analysis of 
experimental data is the excellent agreement between the two approaches as long as 
n\ 2 < 2. This agreement holds for the temperature regime {k^T < 2huj z ) and interaction 
strength (g < 0.15) relevant for current experiments. The Hartree-Fock approach is 
therefore well suited for fitting the wings of the experimental density profiles of a quasi- 
2d gas to extract the temperature and chemical potential. 



7. Probing 2d atomic gases 

This section is devoted to the presentation of some methods that have been used 
for the experimental study of 2d Bose gases. We start with the conceptually simplest 
approach, which consists in the measurement of the steady-state distribution of atoms 
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in a trap. We then turn to the information that can be acquired in a Time-of-Flight 
expansion. Finally we discuss two schemes that give access to the phase coherence of the 
gas. 

7'1. In situ density distribution. - Conceptually the simplest information that can 
be obtained on a 2d gas is a picture of the sample along the direction that is strongly 
confined. Since this degree of freedom is supposed to be frozen out, there is no loss of 
information due to integration along the line-of-sight. This is in contrast to what happens 
in 3d, where one has to resort to a non trivial transformation to reconstruct the spatial 
distribution [98] (see also [99] for a review). 

Assuming that local density approximation (LDA) is valid, the density distribution 
in the trap n(r) can be obtained from the equation of state of the homogeneous system. 
The general form of this equation of state is 

(84) n\ 2 = F(p, k B T, 02) , 

where F is at this stage an unknown function and a-i is the 2d scattering length. Within 
LDA the density n(r) in the trap is calculated by replacing fi by fi — V(r), where V(r) 
is the trapping potential. 

In the quasi-2d regime that is of practical interest (a s <C a z ), we have seen in subsec- 
tion 6'2 that the interactions in the gas are characterized to a good approximation by 
the dimensionless number g = ^/8na s /a z <C 1. In this case eq. (84) can be simplified 
using dimensional analysis; the expression of the phase space density D = n\ 2 must take 
the functional form 

(85) D = G(a,g) with a = . 

Kb 2 

For a gas that is trapped in a harmonic potential mcj 2 r 2 /2, the in situ density profile is 
then given by 

(86) n{r) x* = G (a-^g) , 

where we set as above mio 2 r\, = k^T. This expression clearly shows a scale invariance for 
a given interaction strength g. Suppose that different density profiles n(r) are recorded for 
various temperatures T and various atom numbers N (hence different chemical potentials 
ji). According to cq. (86) the profiles can all be superimposed on the same curve G(a, g), 
provided they are plotted as a function of r 2 jr 2 -, and translated along the x axis by the 
dimensionless quantity a = n/k&T. This scale invariance behaviour has been checked 
with excellent accuracy by M. Holzmann and W. Krauth using quantum Monte-Carlo 
simulations [100]. These simulations were performed for g = 0.15, which corresponds to 
the interacting strength in ENS experiments with Rb atoms. 

For ^ <C 1, various asymptotic forms of the function G(a,g) have been given earlier. 
When interactions can be neglected (g = 0), the equation of state is (cf. eq. (6)): 
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D = — ln(l — e a ). In presence of interactions and for small phase space densities, the 
mean-field Hartree-Fock method amounts to replace fi by n — 2gn into the ideal gas 
result, which leads to the implicit equation D = — ln(l — e Q ~^ D//71 "), from which one 
can extract D as a function of a and g (see also eq. (70)). In the strongly degenerate 
limit, where fi ^> k^T and D 1, density fluctuations are strongly reduced and one 
expects /i = gn, which can be written as D — {2-K/g)a. In the intermediate regime, 
in particular close to the BKT transition point, one can use the results of the classical 
field Monte-Carlo analysis of [44]. The resulting function D = F(a,g) is represented in 
fig. 3 for g — 0.15. The two asymptotic regimes that we just described are indicated 
with dotted and dash-dotted lines. A remarkable characteristic of the function F{a 1 g) is 
precisely the absence of significant features at the critical point for the BKT transition 
(corresponding toaw 0.2 and D « 8 for g — 0.15). This is due to the infinite order of 
the BKT transition, that does not cause any singularity in the dependance of the total 
density n on T or fi. On the other hand the superfluid density n s (plotted as a dashed 
line in fig. 3) is discontinuous at the transition point, but this quantity is not directly 
accessible from an in situ measurement. For a detailed comparison between the results of 
the mean-field Hartree-Fock approach and those obtained from a quantum Monte-Carlo 
simulation and from a renormalization group treatment, see [96] and [63], respectively. 

Finally we note that the analysis of individual images requires a proper knowledge of 
the temperature and the chemical potential. These are usually obtained by fitting the 
wings of the distribution with the appropriate function for the quasi non-degenerate gas 
(see the discussion following eq. (83)). An interesting alternative consists in using in situ 
density fluctuations to determine these thermodynamic quantities [101]. This promising 
method that relies on the fluctuation-dissipation theorem for a non-uniform system has 
not yet been implemented experimentally for a 2d Bose gas. 

7'2. Two-dimensional Time- of- Flight expansion. - Generally speaking, a Time-of- 
Flight (TOF) procedure consists in switching off abruptly the potential confining the 
atoms, letting the cloud expand for an adjustable time, and then measuring the density 
profile. If the role of interactions is negligible during the expansion, the density profile 
after a long TOF is proportional to the in-trap momentum distribution. For a two- 
dimensional system, two types of TOF can be considered. One can switch off the potential 
confining the atoms in the xy plane, while keeping the strong confinement along the frozen 
direction z; we will call this procedure a "2d TOF". Alternatively, one can switch off 
simultaneously the potential in the xy plane and the confinement along z, corresponding 
to a "3d TOF" . We discuss 2d TOF in this section, and 3d TOF in the following one. 

We consider here the case of an isotropic harmonic trap in the xy plane V(r) = 
mui 2 r 2 /2. A 2d gas is initially at thermal equilibrium in this trap, with a density profile 
n eq (r). Suppose that this potential is suddenly switched off at time t = whereas the 
confinement along the z direction remains unchanged. Using the Bogoliubov approach, 
it was predicted in [102] that the subsequent evolution of the density distribution is given 
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Fig. 3. - (a): Phase space density as a function of a = fi/ksT for g = 0.15. Continuous line: 
Total phase space density D = nA 2 ; dashed line: superfluid phase space density n s \ 2 . The 
dotted and dash-dotted lines represent the asymptotic regimes for low and high phase space 
densities, respectively. (b)In situ density profiles in a trap deduced from the left panel using 
the local density approximation. The plot is made for fi/ksT = 0.5 so that the Thomas-Fermi 
radius rrF = \/2/i/maj 2 is equal to rr- 



by the scaling law: 



This means that the global form of the spatial distribution is preserved during the TOF. 
As the expansion proceeds, the interaction energy that was initially present in the gas is 
converted into kinetic energy in such a way that the density profile at time t is obtained 
using a scaling transform of the initial one. We emphasize that this remarkable result is 
stronger than its 3d counterpart [103, 102] which holds only in the Thomas-Fermi regime: 
In the 2d case the scaling behavior is valid both for the superfluid component and for 
the thermal Bogoliubov excitations. This scaling behavior has been recently observed by 
the ENS group [104]. 

The scaling invariance in the expansion of a 2d interacting gas has been explained 
in terms of the SO(2,l) symmetry group for a whole class of interaction potentials U(r) 
between two particles [105]: it is sufficient that U is a homogeneous function of degree 
2, U(ar) = U(r)/a 2 . When this is the case, the result of eq. (87) holds for an arbitrary 
initial state of the 2d gas, irrespective of its temperature. The 2d contact interaction 
potential, which is implicitly assumed in eq. (16), belongs to this class of functions. We 
note however that a true contact interaction is singular in 2d because it leads to ultra- 
violet divergences at the level of quantum field equations. A real interatomic potential 
has a finite range which provides a UV cut-off that eliminates the divergences. This 
regularization will occur if one uses the more precise treatment of atomic interactions 
given in cq. (75). It will lead to deviations with respect to the universal law (87), which 
remain to be evaluated and characterized. 



(87) 



n ( r , t) = Vt n cq {r]tr) , r] t = (1 + uj 2 t 2 ) 1 
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7'3. Three-dimensional Time- of- Flight. - In a 3d TOF both the trapping potential in 
the xy plane and the strong confinement along the z direction are switched off simulta- 
neously. The physics is then very different from that of a 2d TOF. Along the initially 
strongly confined direction z, the atom cloud expands very fast since the momentum 
width Ap z <~ h/Az is large. If the atoms are initially in the ground state of a harmonic 
potential with frequency oj z along this axis, the extension of the cloud is multiplied by 
\/2 in a time t — lu^ 1 . The time scale for the expansion of the gas in the xy plane is 
much longer; it is given by ui~ , where ui <C cu z is the trapping frequency in this plane. 
Therefore it is a good approximation to decompose a 3d TOF into two phases. During 
the first phase, whose duration is a few uj^ 1 (typically 1 ms if cj z /2tt — 3 kHz), the 
thickness of the gas along z increases by a factor much larger than 1, but the xy spa- 
tial distribution is nearly not modified. At the end of this first phase, the interactions 
between atoms have become negligible. During the subsequent phase the expansion in 
the xy plane becomes significant, but on a much longer time scale. It corresponds to the 
expansion of an ideal gas, whose initial state is equal to the state of the system in the 
xy plane before the beginning of the TOF. 

We now focus on the evolution of the xy degrees of freedom during the second phase, 
which is essentially governed by single particle physics. The evolution of the density 
distribution in the xy plane can be determined from the initial one-body density matrix 
gi(r,r') — (r\p^\r'), or from its Fourier transform U(p) with respect to the variable 
r — r' , which represents the momentum distribution in the xy plane. 

In absence of any extended coherence in the gas, gi{r, r') tends to zero when \r — r'\ 
increases, with a characteristic decay length given by the thermal wavelength A. The 
corresponding momentum width is Ap ~ h/X and the spatial distribution after TOF will 
reflect the initial momentum distribution if the TOF duration t is such that Ap t/m 3> 
tt, where tt is the initial size of the gas. For a harmonic confinement in the xy plane, 
this "far field" regime corresponds to ivt 3> 1. Taking u/2tt — 30 Hz as a typical value, 
the far field regime (say u)t > 3) is reached for t > 15 ms. This corresponds to a typical 
value for TOF experiments, which thus give access to the momentum distribution in this 
non (strongly) degenerate regime. 

The situation is very different if a significant condensed fraction is present in the gas, 
as expected in the vicinity and below the BKT transition temperature. In this case we 
have seen in § 5'4 that the size r c of the coherent region of the cloud is r c ~ ry. The 
momentum width Ap c = h/r c of this coherent component is then extremely narrow, 
and it would require a very long TOF to reach the 'far field' regime for this coherent 
component. Taking r c — tt as a typical value, we find that the time t required for a 
significant expansion of this component, i.e. Ap c t/m = r c , is such that tot — k B T/hu. 
For w/27r = 30 Hz and T = 100 nK, this gives t > 300 ms, which is too long in practice 
for a TOF. 

Therefore in the regime where a relatively strong coherence of the gas is present, a 
3d TOF of a realistic duration gives access to a hybrid information. The high energy 
fraction of the gas is in the far field regime and the wings of the density profile after TOF 
give access to the large momentum part of the initial state. On the contrary the central 
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feature corresponding to the condensed, supcrfluid fraction, has not yet undergone a 
significant expansion. The detailed study of the border between these two components 
is still a matter of debate. In experiments with rubidium atoms [27, 106], the density 
profile after a 3d TOF is well modeled by a two-component distribution and fits with 
the line of reasoning we just presented. In contrast, in the experiments performed at 
NIST with a sodium gas [28], an intermediate third component was introduced in order 
to obtain a good description of the density profiles after a 3d TOF. This component 
corresponds to a phase with a spatial scale of coherence that is intermediate between the 
microscopic length A and the macroscopic one r T , and it is qualified as a "non-supefluid 
quasi-condensate" in [28]. 

It is interesting to note that 3d TOF is the most common and natural experimental 
method used in the studies of 3d atomic gases. However, in hindsight, its availability is a 
non-trivial and rather serendipitous feature of atomic systems for studies of BKT physics. 
In combination with the finite-size induced condensation, the ability to suddenly turn off 
the interactions through the fast z-expansion provides a much more striking signature of 
the BKT transition [27, 28, 106] than one might have theoretically expected. Thinking 
strictly in 2d, the transition is extremely smooth and one would not naturally expect 
to see such a dramatic signature in any quantity except the superfluid density. As we 
discussed in § 7'2, in 2d TOF the observed density distribution indeed varies smoothly 
across the transition. 

So far we have discussed the 'average' density profile in 3d TOF, which theoretically 
corresponds to the average of a large number of images obtained under same conditions. It 
is also interesting to consider the density noise in individual images, which can be related 
to the phase noise of the gas before expansion. This connection has been exploited for 
quasi-ld gases since 2001 [107]. For the 2d case, it has been shown theoretically in [108] 
that the two-point density correlation function after TOF can provide information on 
the in situ gi function, at least in the superfluid regime. This method is also specific 
to 3d TOF, where the phase noise evolves into density noise in interaction-free ballistic 
expansion. 

7'4. Interference between independent planes. - Since an important aspect of the 
physics of 2d Bose gases is related to phase properties, it is natural to investigate mea- 
surement schemes based on interferometry. We start with the proposal by Polkovnikov 
et al. [109] which showed how a single experimental procedure could characterize both 
the normal regime (exponential decay of gi) and the supcrfluid regime (algebraic decay 
of g\) (see also [110] for a more complete review). Consider two independent, infinite 
planar gases located at z a = +d z /2 and z b = —d z /2. They are prepared in identical 
conditions, i.e. they have the same temperature and the same density. We perform a 3d 
time-of-flight of duration t, that is chosen such that the final extension along z of each 
cloud is large compared to the initial separation d z between the planes. The two clouds 
thus overlap and we want to extract information about the one-body correlation function 
gi from their interference pattern (fig. 4a). 

The state of each plane is described by the wave function ip a /b{x, y). After expansion, 
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Fig. 4. - (a) Principle of an experiment giving access to the interference between two independent 
planar gases, observed after time-of- flight, (b-d): Examples of interference patterns measured 
with the experimental setup described in [26]. The imaging beam is propagating along the y 
axis. The pattern (b) is obtained with very cold gases, whereas (c) corresponds to a larger 
temperature. The dislocation in (d) is the signature for the presence of a vortex in one of the 
two gases. 



the spatial atomic density n is modulated along any line parallel to the z axis with the 
period D z = ht/md z [111]: 



For simplicity we have omitted in the above equation a global envelope factor giving 
the variation of the density along the z axis. Also we have neglected the expansion in 
the xy plane during the TOF. As explained above there exists a range of TOF duration 
where this is valid, if the trapping frequency uj in this plane is much smaller than uj z . We 
see from eq. (88) that the local (complex) contrast of the density modulation is ip a ipt- 
Experimentally one cannot measure this quantity along a single line, and one rather has 
access to the average contrast over a region of finite area A in the xy plane. In particular 
if one performs absorption imaging along the y axis (fig. 4b-d), the image involves an 
integration of the local contrast ip a ipl along the y direction ( 8 ). Averaging the result of 
this contrast measurement over a large number of realizations, one can define the average 
contrast C(A) : 



(89) 



C 2 {A) 



1 

A?' 



ip a (r)ipt(r) d 2 r 



Using the fact that the fluctuations of the wave functions ip a and ipb are uncorrelated 
and taking advantage of the translational symmetry of the system, we find: 



(90) 



c2 ( A ) = \J A \9i(r)\ 2 d 2 r , 



( s ) The length over which the line-of-sight integration occurs can be adjusted by a proper 
"slicing" of the cloud just before the imaging process, as in [111]. 
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where gi(r) = (ip*(r)tpj(0)} for j = a,b. Suppose for simplicity that the area A is a 
square and consider the two cases of an exponentially decaying gi(r) oc e~ r / e , with a 
characteristic length I <C \[A (normal fluid), and an algebraically decaying g\(r) oc r~ Q , 
with an exponent a < 1/4 (supcrfluid regime). In the first case, the integral is nearly 
independent of A and C 2 (A) oc A^ 1 . In the second case we find C 2 (A) oc A~ 2a which 
corresponds to a decay always slower than A~ x / 2 . This method is very appealing in the 
sense that the measurement of a single number, i.e. the exponent r\ characterizing the 
variation of C oc A~ v , is sufficient to identify the two possible regimes of a 2d Bose gas, 
and obtain the value of n s X 2 — 1/r] in the superfluid case. 

A measurement scheme inspired by this method was implemented at ENS [26] , and it 
indeed revealed a relatively rapid variation of the exponent 77, in qualitative agreement 
with what is expected near the BKT transition point in the center of the trap. However 
some notable deviations with respect to the original proposal must be stressed. First, 
the measurement was performed with anisotropic samples, with lengths L y <C L x . The 
imaging beam was propagating along y and the measured contrast involved a line-of-sight 
integration over the full length L y ( 9 ), which formally breaks the translational invariance 
that we used to prove eq. (90) . Also the presence of a trapping potential in the experiment 
causes an additional softening of the transition, because of the inhomogencity of the 
density along the line-of-sight of the imaging beam. Finally we note that even deep in 
the superfluid regime where n s X 2 >• 1, the anisotropy of the sample adds some complexity 
as discussed in §4'5. At large distances (A^ > n s \ 2 L y ), g\ starts to decay exponentially, 
which complicates the analysis of the dependence of C 2 on A x . In summary the rapid 
increase in coherence that occurs in the vicinity of the BKT transition point is sufficiently 
robust to be revealed experimentally in the average contrast of the interference pattern, 
but it is difficult to provide a quantitative analysis of the experimental measurements for 
the variations of C 2 (A) over a large range of A x . 

A subsequent experiment at ENS has compared the conditions for observing a signif- 
icant interference contrast between the planes and for measuring a clear bimodal density 
profile after a 3d TOF [27]. The onsets of the two phenomena were found to coincide 
within experimental error. Furthermore the spatial part of the gas that gives rise to a 
visible interference signal coincides with the central, "non expanding" component of the 
TOF profile. 

An important outcome of the experiments on the interference between two planes is 
a direct evidence for thermally activated vortices. At low temperatures, long-wavelength 
phase fluctuations (phonons) result in smooth variations of the phase of the interference 
fringes, such as seen in fig. 4c. However, if a single isolated vortex is present in one of the 
two planes while the phase profile of the other plane is smooth, the interference pattern 
exhibits a sharp dislocation at the coordinate x of the vortex core. Such dislocations 



( 9 ) The area A was varied by changing the integration distance along x. In this case, the 
BKT transition causes a crossover from r\ = 1/2 for an exponentially decaying g\ function (with 
a decay length I <§C A x ), to 77 = 1/4 for a superfluid state. 
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have been observed experimentally [81, 26, 69] and an example is shown in fig. 4d. 
The occurrence probability of these dislocations has been measured as a function of 
temperature [26]. The number of dislocations increases with T, until one reaches the 
temperature at which no interference is visible anymore. Moreover, the relatively sharp 
increase in the probability of dislocations experimentally coincides with the increase in 
the exponent r\ characterizing the decay of g\ [26]. Such dislocations also appear in 
a classical field simulation mimicking the interference between two planar gases [40]. 
They result from the thermal activation of a vortex pair for which the two members 
are sufficiently separated from each other. In principle one should also observe in the 
interference patterns tightly-bound vortex pairs where the two members are separated 
by <~ £. However these pairs only shift the fringe pattern by a small fraction of the fringe 
spacing, which is below the current sensitivity of the experiments. 

7'5. Interfering a single plane with itself . - An interesting alternative to the two-plane 
interference described above consists in preparing a single plane of atoms and looking at 
its "self-interference" , using a Ramsey-like method [28] . The gas is initially prepared in 
an internal state |1). Half of the atoms are coherently transferred into another internal 
state 1 2) by a stimulated laser Raman process (tt/2 pulse) that also provides a momentum 
kick k to the atoms. After this process, the part of the cloud in |1) is still globally at 
rest and the part in |2) moves with the global velocity v = hk /m. After an adjustable 
time t a second tt/2 Raman pulse remixes the amplitudes of |1) and |2) and provides 
a momentum kick fco — fei. Immediately after this second Raman pulse, one measures 
the spatial density distribution in |2). This distribution exhibits a modulation along the 
direction ki, resulting from the interference between the initial state of the cloud and 
the state displaced by the distance R = v t: 



Note that we assume here that no collision occurred during the time t between the part 
of the cloud at rest in |1) and the part moving at velocity v in state |2). This is a 
valid assumption for the weakly interacting sodium gas of [28]. The modulated density 
profile in eq. (91) gives a direct access to the function g\(r, r — R). It can be observed 
with an imaging beam along the z direction, so that its measurement does not involve 
any line-of-sight integration. This method can then reveal finer details than the one 
presented in § 7'4. In particular the authors of [28] could observe a gradual increase of 
the coherence length I of the cloud, as expected from eq. (56). For small phase space 
densities the measurement gives I <~ A, and I increases to much larger values when 
the temperature decreases towards the critical temperature Tbkt- When T < Tbkt 
a significant interference contrast is observed for all values of R within the size of the 
central superfluid region. 



(91) 



n(r) cx \ip(r)\ 2 + \ijj(r - R)\ 2 + (ip(r)tp*(r - R)e 



+ C.C.) . 
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8. Conclusions and outlook 

We have reviewed in these notes the theoretical basis for the understanding of the 
physics of 2d quantum fluids, and discussed some recent experiments performed with 
atomic gases. These experiments have given access to some aspects of 2d physics that 
were previously hidden or not measurable in other physical systems, such as the exis- 
tence of thermally activated individual vortices or the spatial variation of the first order 
correlation function gi(r). However a number of issues is still open in the physics of 2d 
quantum gases, and we outline below some topics which are likely to be of experimental 
and theoretical interest in the future. 

Higher-order correlation functions. The matter-wave interference between two statisti- 
cally similar, but independent quasi-condcnsates (such as shown in figure 4), can reveal 
a wealth of information on the correlations within each individual 2d gas. So far only a 
fraction of this information has been harnessed, with the study of the average contrast 
of the interference pattern integrated over some area of interest. A convenient tool for 
extracting more complete information on g\ as well as higher-order correlation functions 
is the full statistical distribution of interference contrasts. Two limiting cases can easily 
be characterized: (i) If the two independent fluids arc fully condensed, each image shows 
a 100% contrast, with the position of the fringes fluctuating randomly from shot to shot, 
(ii) If each cloud exhibits only short-ranged correlations, the observed interference results 
from many uncorrelated fringe patterns along the light of sight, and the distribution of 
contrasts is an exponential function. For Id gases, it is possible to describe quantitatively 
the transition between these two limiting cases [112], and the experimental results [113] 
are in good agreement with the predictions. In the 2d case, the evolution of the contrast 
distribution through the BKT transition is still an open problem. 

Out- of- equilibrium dynamical effects. Throughout this paper we restricted our discussion 
to the equilibrium properties of a 2d Bose fluid. The study of dynamical effects, such 
as transient regimes, can reveal additional information about the system. For example 
Burkov et al. [114] have studied the dynamics of decoherence between two planar Bose 
gases, assuming that their local phases are initially locked together, and then the two 
gases are allowed to evolve independently This can be achieved experimentally by having 
a weak potential barrier and hence large tunnel coupling between the two planes for 
t < 0, and then suddenly raising the barrier at t = 0. The contrast of the interference 
between the two gases gives access to the evolution of the phase distribution under the 
influence of thermal fluctuations. In [114] this contrast was shown to decay algebraically 
at long time, as i - ^, with the exponent ( proportional to the ratio T/Tbkt- Therefore, 
in addition to being a stringent test of thermal decoherence in a quantum many-body 
system, this out-of-equilibrium study could constitute a novel thermometry method. A 
related phenomenon occurs in Id systems, where the interference contrast is predicted to 
decay as exp(— (t/t a ) 2 ^ 3 ) (with t constant) [114], and this prediction is nicely confirmed 
in the experiments by the Vienna group [115]. 
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Transition from 2d to 3d behavior. The possibility to vary the tunnel coupling between 
two or more planar gases can also be used to study the so-called "deconfinement tran- 
sition" [116], corresponding to a gradual evolution from 2d to 3d behavior. The phase 
coherence between the planes will build up as the strength of the coupling is increased, 
creating a situation that is reminiscent of the high-T c cuprate superconductors. For 
a large number of parallel planes, the deconfinement transition should give rise to a 
true Bose-Einstein condensate [116]. The two-plane situation is also very interesting, 
and can lead to the observation of the Kibble-Zurek mechanism [117]: the superfluid 
transition temperature is higher for two coupled planes than for a single one, so that 
sudden switching on of the coupling between the planes (initially in the normal state 
but close to the single-plane critical temperature) constitutes a quench of the system, 
and one could observe the subsequent dynamical apparition of a macroscopic quantum 
(quasi- ) coherence . 

Tunable interactions. As we have seen throughout the paper, interactions between par- 
ticles play a crucial role in our understanding of the superfluid phase transition and 
condensation in 2d fluids. In contrast to the conventional 3d BEC of an atomic gas, 
where the critical temperature can to a good approximation be predicted using the ideal 
gas model, the BKT transition is fundamentally interaction-driven. The strength of in- 
teractions also affects a variety of other phenomena such as the suppression of density 
fluctuations in the normal state and the connection between the 2d Bose fluid and the 
XY model. It would therefore be interesting to revisit the various effects described in 
these notes while continuously tuning the strength of interactions with a Feshbach reso- 
nance [48, 49]. In the weak coupling regime (g < 10" 1 ) we expect a gradual change from 
the BKT-dominated to the BEC-dominated behavior, as discussed in Sections 4 and 5. 
Further, it would be very interesting to explore the strong-coupling regime (g > 1), which 
is closer to liquid helium films. This regime, which is outside the domain of validity of 
the Monte-Carlo results [29, 44], corresponds to the case where the scattering length a s 
becomes comparable the thickness of the sample along the kincmatically frozen direction, 
a z (see 6'2). There the very nature of two-body interactions is expected to change from 
3d to 2d [84, 85, 118, 119]. Therefore, experimentally reaching the condition a s > a z 
would correspond to producing a "truly 2d" as opposed to a quasi-2d Bose gas. 

Superfluid density. Generally speaking, studies of coherence and correlation functions in 
a 2d fluid, which are well suited to experimental tools of atomic physics, are a natural 
complement to the "traditional" studies of superfluidity based on transport measure- 
ments, which are well suited to other physical systems such as liquid helium films [2]. 
For example, we have so far assumed that the two types of measurements probe the same 
superfluid density (see e.g. § 7'4). However this correspondence may in fact depend on 
the theoretical model and the exact definition of the superfluid density, and be valid only 
within the effective low-energy theories. It is therefore important to stress that super- 
fluidity in the traditional transport sense has not yet been directly observed in atomic 
2d Bose gases (see e.g. [116]). Establishing atomic 2d gases as experimental systems in 
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which both coherence and transport measurements of superfluidity could be performed 
would be an important advance, as it would allow experimental scrutiny of the theoreti- 
cal connections between the two types of probes, and a direct comparison of the different 
definitions of superfluidity. Two promising schemes for a direct measurement of the su- 
pcrfluid density (as traditionally defined through transport properties [46]) in an atomic 
gas have recently been proposed [120, 121]. The first scheme [120] is based on extracting 
the superfluid density from the in situ density profiles of a rotating 2d gas. The second 
scheme [121] is based on using a vector potential generated by Raman laser beams to 
simulate slow rotation of a gas [122], and allows direct spectroscopic measurement of the 
superfluid density. 

Note added in proof . A measurement of the equation of state of a 2d Bose gas for various 
interaction strengths has just been reported in [123]. 
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